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ABSTRACT 


A series of low-speed wind-tunnel investigations were conducted to determine 
the aerodynamic behavior of a grooved inclined cylinder representing a long trailing 
wire antenna towed from an orbiting airplane. The large angle-of-attack range of the 
trailing wire required two different model configurations. The first configuration, 
using full-scale wire lengths suspended between steel stanchions, was mounted on a 
flush four-degree-of-freedom wall balance. The second configuration used a 15-scale 
grooved cylinder model with an ogive nose mounted on a six-degree-of-freedom sting 
balance. Wall balance wire data, valid for higher angles of attack, were integrated 
with low angle-of-attack sting balance data. Empirical relationships for the normal 
and axial force coefficients were verified with historical references for tested clean 
circular cylinders and extended for the grooved configurations. Existence of a side 
force coefficient due to circulation caused by the helical grooves was discovered, 
expressed analytically, and verified with flow-visualization techniques. Finally, the 
experimental coefficients were used to improve an existing simulation model 
describing the static equilibrium conditions of a cable towed by an airplane in a 
circular orbit. Inclusion of the side force influence in the static model proved 


consistent with the lateral skew angle and direction observed during flight test. 
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I. INTRODUCTION 


The use of wires in the field of aeronautics started well before the Wright Flyer 
took to the air. Initially, wires were used as structural support for bi-wing gliders, 
and their use extended to powered biplanes. With their use came experimental 
Studies and tests, and the flow over an infinite cylinder became the baseline for many 
disciplines in classical aerodynamics and fluid mechanics. 

With the rapid advance of aircraft design, and the development of the 
monoplane with wing spars, metal skins, and monocoque fuselages, the use of wires 
and cables became restricted to internal flight control systems, and a new purpose, 
as instruments to tow other aircraft. Interest in towing cable form and stability 
manifested as early as 1934 with a study by Glauert [Ref. 1]. 

With the advent of radio and the use of a long trailing wire as an antenna, 
concern for wire stability and form increased. The Navy EC-130, circa 1970, was the 
first of two airplanes designed to trail 15,000 to 25,000 feet of wire in a circular orbit 
for use as a very low frequency (VLF) antenna. The wire, with a cone-shaped 
drogue at the end, exhibited extreme oscillations in shape and tension. These 
oscillatory problems continued with the 1987 replacement airplane, the E-6A. The 
trailing wire came in two configurations: a 19-strand wire wrapped by a flattened 
copper sheath (1x19), and three 7-strand bundles twisted in a triple thread helical 


groove of 15° (3x7). 


Three major discrepancies were documented during E-6A flight test of the 3x7 
wire. The first discrepancy from the oscillatory behavior was contact between the 
wire and towplane’s horizontal tail. The conical motion of the first 50 feet of wire 
and subsequent airplane contact has potential to entangle horizontal flight control 
surfaces. The second discrepancy stemmed from wire failure due to oscillatory 
tension forces exceeding the ultimate strength of the wire. The third discrepancy was 
the mission degradation displayed by the reduced transmissivity of the wire during 
oscillations. Also discovered during flight test was an exhibited trait of the 3x7 wire 
to trail approximately 10° to 15° right of centerline (looking aft) while in level flight 
and zero degrees bank angle. 

Limited simulation studies and expensive flight test experimentation were 
conducted over the past 20 years to model and control trailing wire oscillations. 
Recent studies by Clifton [Ref. 2] currently provide the best system model and viable 
solutions to the nonlinear oscillatory dynamics presented by the 3x7 wire. Clifton’s 
model commenced with the establishment of no-wind static equilibrium conditions 
of the 3x7 wire, towplane, drogue system while in a circular orbit. The simulation 
model was extended to include dynamic conditions and a wind forcing function. 
Accuracy of the Clifton model was only limited by the use of approximate force 


coefficients for the 3x7 trailing wire and drogue. 


The purpose of this research was to: 


@ experimentally determine all aerodynamic traits of the E-6A trailing wire and 
drogue configurations, 


@ discover, determine, and explain the existence of a side force induced from the 
helical grooves of the 3x7 wire at low angles of attack, and 


@ establish the amplitude of the side force, and by refining the Clifton static 
simulation model, determine its effect on the E-6A antenna system. 


Il. BACKGROUND 


A. BASIC AERODYNAMIC THEORY 

The study of the infinite cylinder has become the foundation of aerodynamic 
analysis. From classical fluid mechanics, stream functions of two elementary flows, 
a doublet and a uniform flow field, can be superimposed to produce a streamline 
describing flow about the surface of a circular cylinder. Pressure fields can then be 
calculated for two-dimensional, invisid flow with analytical solutions determined. 
Integration of the pressure field in an invisid, irrotational flow about a cylinder 
results in forces in the direction of flow, drag, and orthogonal to the flow, lift, equal 
to zero. By adding circulation to the flow fields, the force orthogonal to the flow 


(lift) acting per unit length was nonzero and expressed by Kutta-Joukowski as 


F = p*(V_xT) (2.1) 
Where: F = force vector 
Po = freestream density 
V., = freestream velocity vector 
I = circulation vector 


A circular cylinder in viscous flow acts as a bluff body producing drag primarily due 
to flow separation over the downstream part of the cylinder. The drag about a 
cylinder, due to separation-induced pressure drag, results in a drag coefficient 


dependent on the nondimensional Reynolds number, Re, defined as 


p_*V_+«d 


Re = ————_ (2.2) 
ar 
Where: P. = freestream density (slugs/ft’) 
V.. = freestream velocity (ft/s) 
d = diameter of cylinder (ft) 
Lo. = freestream viscosity (slug/ft*s) 


The relationship between C, of a circular cylinder and Reynolds number is presented 


in Figure 2.1. [Ref. 3] 
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Figure 2.1 Circular Cylinder C, vs Reynolds Number 


In the subcritical Reynolds number range, from 4x10° to 3x10’, the C, of the cylinder 
is an approximately constant value of 1.2. Subcritical flow over a cylinder is 
characterized by laminar conditions until separation around the +90° from 
Stagnation. The wake is characterized by an asymmetric "vortex street" first analyzed 
by von K4rm4n. The flow reaches critical Reynolds number around 3x10° as 


transition from laminar to turbulent conditions occur. Supercritical flow is 


characterized by a lower value of Cp due to turbulent layer separation past + 100° 
producing a smaller wake and less form drag. [Ref. 4] 

An experimental study of smooth and stranded wires was conducted as early as 
1917 by Relf and Powell [Ref. 5]. In 1951 Hoerner published Aerodynamic Drag 
followed in 1958 by Fluid Dynamic Drag. His works presented empirical 
relationships for the lift and drag coefficient of inclined wires, cables and cylinders 
at subcritical Reynolds numbers. In a body axis system for an inclined wire as shown 
in Figure 2.2, the normal force coefficient, Cy, is determined from the force 
perpendicular to the wire, and the axial force coefficient, C,, from the force parallel 
to the wire. Stability axis lift and drag coefficients can be resolved by transforming 


Cy and C, using 


A 


C ~j C 
L _ COSa@ Sina : N (2.3) 
on sina cosa | {C 
Hoerner produced an analytic expression from empirical results based on Cy, only. 


Summarized: 


o aaa 
Cy = Co pasic *SIN" & 





C, = Cp pasie *Sin’ a *cosee (2.4) 
C 
C pase. oe = 
; sin* a 


N/(Q*S) 

normal force (lb,) 

dynamic pressure = 0.5*p,,*V,.” (Ib,/ft’) 
reference area (ft’) 


Where: 


nO Zz! 
lh mom 


Setting Cp pasic = 1.1, Hoerner’s closed-form solutions are presented in Figure 2.3. 
Hoerner did not consider an axial force but added a frictional term AC, = rC,=0.02. 


[Ref. 6] 
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Figure 2.2 Inclined Wire Coordinate System 
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Figure 2.3 Hoerner’s Empirical Relationships [Ref. 6] 


In 1970, Bootle [Ref. 7] extended Hoerner’s expressions to supercritical 
Reynolds numbers at low Mach numbers. Further studies concentrated on critical 
and supercritical flow over slender, axisymmetric, finite-length bodies in response to 
interest in missile related aerodynamics. Studies by Alnosnino and Rom using a sting 
balance cone-cylinder model discussed the magnitude of side force coefficient Cy at 
higher angles of attack, and the effect of symmetric blowing on alleviation of side 
force (Ref. 8]. Vortex-induced side force investigations continued with high-angle-of- 
attack studies by Reding and Ericsson [Ref. 9]. Their experiments postulated that 
the maximum side force occurred at critical Reynolds numbers on an order of 
magnitude equal to the ogive cylinder normal force. Other studies have addressed 


vortex separation points by the measurement of skin friction coefficient (Ref. 10]. 


Currently, research only addresses vortex induced separation and side forces, as 
applied to missiles, at high angles of attack and high Reynolds numbers. 

This paper addresses the side forces generated by helical grooves on an infinite 
cylinder, representative of a: wrapped inclined wire, at low angles of attack and 


subcritical Reynolds numbers. 


B. 3x7 WIRE 

One towed cable used for long trailing antenna applications is the 3x7 wire. 
Constructed of three sets of six symmetrically placed steel wires around a seventh 
Strand, the 0.158-inch-diameter wire was wrapped in a triple helix with a 1.87-inch 
pitch. The forming process of the three wire groupings yielded a helical indentation 
Or groove with an approximate helix angle of 15°. Shown in Figure 2.4 and 2.5, the 
3x7 wire exhibited failure at a 3000 |b, load. With a length up to 25,000 feet, the 3x7 
wire flew in an aerodynamic environment with altitudes from 3000 to 20,000 feet, 
airspeeds from 36 knots (60 ft/s) to more than 225 knots (375 ft/s), and angles of 


attack ranging from 10° to 90°. 
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Figure 2.4 3x7 Wire 





Figure 2.5 3x7 Wire (End View) 
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C. DROGUE 

The antenna system drogue was designed to provide aerodynamic stability and 
adequate tension during reel-in and reel-out operations. The 31.7-inch-long drogue 
has a base diameter of 24 inches and a weight of 87 Ib,. The center of gravity was 
measured 13.0 inches aft of the 22.5° apex. In actual operation, the wire is attached 


to the forward end of the hollow-cone, weighted-nose drogue. 


D. STATIC CONDITIONS SIMULATION COMPUTER MODEL 

The static conditions simulation program, TAC17, was the foundation for the 
dynamic, nonlinear system simulation model by Clifton [Ref. 2]. The static model 
simulated an airplane in circular orbit at a specified altitude towing a trailing wire 
in no-wind conditions. Although the program allowed freedom of choice for airplane 
altitude, velocity, and bank angle, as well as wire length, baseline airplane 
operational conditions were 18,325 feet of altitude, 156 KEAS (knots, equivalent 
airspeed), and 34° bank angle towing 20,290 feet of wire. The baseline conditions 
corresponded to a particular set of flight data used by Clifton for correlation 
purposes. Program design for an orbiting airplane required a minimum bank angle 
of 3° (left) to preclude an infinite orbit radius. Wire and drogue aerodynamic 
coefficients were defined values, and the model used central difference numerical 
methods with the wire represented by 200 equi-length segments. Program operation 
commenced from an initial guess position of the drogue, and calculated up the 200 


grid points to the towplane position. Calculations were iterated until towplane 


ji 


boundary conditions (inputed) were satisfied to a specified error. Outputs from the 
program included radial, theta, and z position, tension, angle of attack, true velocity, 
and Reynolds number for each wire grid point. A detailed discussion of the static 


conditions simulation model genesis and operation is presented in Reference 2. 


U 


IT. EXPERIMENTAL SETUP 


A. WIND TUNNEL 

Experimentation was conducted in the Naval Postgraduate School (NPS) 
horizontal low-speed wind tunnel. Presented in Figure 3.1, the Aerolab single-return, 
closed-circuit tunnel was powered by a 100 hp electric motor driving a three-blade, 
variable-pitch fan. A four-gear transmission and a 10:1 contraction ratio allowed for 
‘test section speeds to 200 miles per hour. Low test section turbulence intensity of 
0.2% was attributed to stator blades immediately aft of the fan, and two fine wire 
mesh screens six inches apart in the settling chamber, in addition to turning vanes in 
each corner. The 8.75 ft’ test section, 45 inches wide by 28 inches in height, was 
slightly divergent to counter the effective contraction caused by boundary layer 
growth. A 5/100 tunnel diameter breather slot, immediately downstream of the test 
section, helped maintain approximate atmospheric static pressure conditions. Swing 
windows on either side of the test section and frosted glass corner fillet fluorescent 
lights provided adequate illumination, visualization, and access to the test models. 

Test section dynamic pressure, Q, was determined by measuring the static 
pressure difference, Ap, between four manifold-flush taps in the test section and a 
similar set of four taps in the settling chamber. Connected via a common manifold, 


the Ap was presented on a micromanometer and digital display. The Ap was 
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Figure 3.1 NPS Horizontal Low Speed Wind Tunnel 


converted into dynamic pressure using the relationship in equation (3.1) calculated 


from a previous tunnel calibration. 


Q = ae - 2.046+K+Ap (3.1) 
Where: dynamic pressure (Ib,/ft’) 

freestream velocity (ft/s) 

Static pressure difference (cm of H,O) 
tunnel calibration constant = 1/0.93 


| 
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Tunnel air temperature was measured using a dial thermometer extending into the 
settling chamber. A more detailed description of the wind tunnel is presented in 


Reference 11. 
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B. BALANCES 


I. Wall Balance 

An external, cylindrical, reflection-plane (wall) balance, shown in Figure 
3.2, was used to measure forces and moments of the wire at the higher angles of 
attack. Flush mounted in the test section, the wall balance was designed to measure 
normal and axial forces to 150 lb, and their respective moments. Built by NPS 
personnel in 1974, the four degree of freedom balance was designed with two pairs 
of strain gage bridges orthogonally mounted on flexure links separated vertically by 
26.5 inches. Each bridge circuit had four active legs for automatic temperature 
compensation. The lower bridge A, and the upper bridge B, each output an axial 
and normal voltage (E,,, E,,, E.,, E,,) which were converted into forces and moments 
using the results of a balance calibration described in Appendix A. Electrical 
schematics of the bridges and pin connections are presented in Reference 12. 

The balance column was rigidly mounted on an electrically-controlled 
turntable capable of rotating -18° to +200° from the centerline orientation shown 
in Figure 3.3. A 15.625-inch-diameter aluminum plate was mounted with eight 
screws atop the reflective plane and flush with the floor of the test section. A 0.125 
inch gap existed between the plate and reflective plane to prevent contact and ensure 
accurate measurements when under load. Models were mounted on the aluminum 


plate for tests. 
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Figure 3.2 NPS Wall Balance 


2. Sting Balance 
A six-degree-of-freedom, one-inch-diameter, Mark XIV, internal Task’ 
balance was used to measure forces and moments for models representing lower 
angles of attack of the wire, and for the drogue. The 4.9-inch-long balance had eight 
3500 wheatstone bridges: two normal force (N1, N2), two side force (S1, S2), two 
axial force (Al, A2), and two rolling moment (11, 12). The two axial force and two 


rolling moment bridges were each connected in parallel to produce 1750 circuits, 
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Figure 3.3 Wall Balance Coordinate System 


and therefore six electrical outputs: N1, N2, S!, $2, A, |. Maximum balance loads 
were 400 Ib, in the normal channels, 200 lb, in the side force channels, 100 Ib, axially, 
and 250 inlb, of rolling moment. On loan from NASA-Ames Research Center 
through the Navy-NASA Joint Institute of Aeronautics, the balance was calibrated 
to a 5 Vy, bridge excitation voltage by NASA-Ames personnel. Additional 
calibration, presented in Appendix B, was required to ensure proper force 
measurement, and moment focal point position prior to testing. Calibration 
constants and their accuracies are presented in Appendix B, Tab I. 

The balance was mounted into a 6.75-inch sleeve extender and locked 


with set screws onto a "U" frame as shown in Figure 3.4. The twenty-four 36-gage 
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Figure 3.4 MK XIV 1" Task (Sting) Balance 


wires were fed through the sleeve and frame, and out the tunnel. Wire slack was 
provided to allow +65° of angle of attack (AOA). Balance orientation is presented 
in Figure 3.5. The balance had four pairs of pin recesses orthogonally positioned for 


model mounting. 
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Figure 3.5 Sting Balance Coordinate System 


C. ACQUISITION SYSTEM 
1. Signal Conditioners / Pacific® Amplifiers 

The electrical bridge outputs from the wall and sting balances were 
connected to individual signal conditioners that controlled bridge excitation. The 
excitation voltages were 10 V,- for the four wall balance channels and a 
precalibrated 5 Vp- for the six sting balance channels. Sting balance channels were 
connected to the cannon plugs on the signal conditioners as shown in Figure 3.6. 
The conditioned signals were amplified with a gain of 1000 by individual Pacific’ 
8255/6 amplifiers and routed to a National Instrument MC-MIO-16L-9, 50 pin 


input/output (I/O) connector presented in Figure 3.7. Wall balance channels E,,, 
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Figure 3.6 Sting Lead/Cannon Plug Connections 


E.. Em Ey, were connected to pins 4/3 (lead/ground), 8/7, 12/11, 16/15, 

respectively. Sting balance channels N1, N2, A, S1, S2, | were connected to pins 4/3, 
6/5, 8/7, 10/9, 12/11, and 14/13 respectively. 
2. Data Sampling / Computer System 

The MC-MIO-16L-9 board has a 12 bit (i.e. 2" bits) analog to digital 

(A/D) converter with a 9 ys conversion time and design acquisition rates up to 100 

Kbytes/s Ref 13}. Digital resolution was 4.88 mV with bridge excitation of 10 Vy. 

(-2048 to 2047 bits) and 2.44 mV with an excitation of 5 V,,-. Data acquisition was 

controlled by programs written and complied in QuickBasic using Microsoft 


QuickBasic 4.5° software and run on an IBM PS/2° microcomputer. National 
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Figure 3.7 MC-MIO-16L-9 I/O Board 


Instruments LabWindows’ interactive software version 1.1 was utilized in the 
programs to command data sampling and averaging. One thousand samples per 
channel were taken for each data point at a rate of 1770 Hz. Time averaged 


sampling required 2.26 seconds for the wall balance and 3.39 seconds for the sting 


balance. 
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The QuickBasic programs converted the averaged voltage samples to 
forces and moments for output and analysis. The wall balance program multiplied 
the four voltage samples by the 4x4 calibration matrix, [K], to produce normal and 


axial forces and moments in balance coordinates. 
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Voltage readings, in addition to the force and moment calculations, were output to 
two files for analysis. The sting balance program utilized nonlinear equations derived 
by Yuan [Ref.14], multiplying calibration constants in Appendix B, Tab I by the six 
voltage readings. Additionally, nonlinear interaction equations were iterated to 
simultaneously solve cross channel dependence of the balance. The primary output 
file consisted of normal, side, and axial forces, and pitching, yawing, and rolling 
moments. Both wall and sting balance acquisition programs are presented in 


Appendix C. 


D. MODELS 


1. Wall Balance Models 
All wall balance models were mounted on the 15.625-inch-diameter 
rotating aluminum plate. The baseline model, shown in Figure 3.8, consisted of six 


(.3125-inch-diameter steel rods horizontally oriented between two 0.625-inch- 
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Figure 3.8 Clean Cylinder (Rod) Model (105° AOA) 





diameter, 13.5-inch-tall steel] stanchions. Each rod presented 12 inches of length to 
tunnel flow. By vertically separating the rod by 2.0 inches, less than 1% interference 
effects from an adjacent rod was calculated. The rods, representing long circular 
cylinders, were used to verify experimental procedures, for analysis, and to provide 
a baseline for the 3x7 wire comparison. 

Six actual 3x7 wires of 0.158-inch diameter were also mounted between 
the 0.625-inch-diameter stanchions. Shown in Figure 3.9 at an angle of attack 90° 
to the flow direction, the wires presented six feet of total length at full scale 
Reynolds number. Finally, two clean 0.625-inch-diameter stanchions were 
constructed for calculation of stanchion forces alone. Sketches of wall balance 


models are presented in Appendix D. 
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Figure 3.9 3x7 Wires (90° AOA) 


Ze Sting Balance Models 

Sting balance models were mounted on the 1.0-inch-diameter Mark XIV 
internal balance. The baseline model, shown in Figure 3.10, consisted of a 21.0-inch 
long, 2.375-inch-diameter aluminum cylinder with a 5.35-inch long, two-caliber 
tangent-ogive nose with radius of 10.094 inches blunted at the nose by an 0.375 inch 
radius. When mounted, the nose declined 1.6° from the horizontal in the (-) side 
force coordinate direction. The balance focal point was 10.0 inches aft of the nose 
cylinder joint. Additionally, 80 and 220 grit sandpaper was mounted at the nose base 
to trip the boundary layer and disrupt nose-generated asymmetric flow separation and 
vortex shedding [Ref. 15]. The clean cylinder sting model was used to verify 
experimental procedures, as a baseline for integration with wall balance data, and 
comparison with helix sting models. 

Helix models were constructed to represent the 3x7 wire. The right helix 


model, shown in Figure 3.11, was a 15-times scale model constructed with the same 
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Figure 3.10 Clean Cylinder Sting Model 


relative 1.87-inch helical pitch and direction of the 3x7 wire. The 21.0-inch-long, 

2.375-inch-diameter grooved section was representative of 1.4 inches of 3x7 wire and 
contained 0.75 revolutions. The left helix sting model, shown in Figure 3.12, was 
identical in groove pitch but opposite in direction to the right model. Nose 
dimensions were identical to the clean baseline cylinder. Finally, a 0.32-scale sting 
model of the wire system drogue is shown in Figure 3.13. Drogue model length was 
9.264 inches with a base diameter of 7.675 inches. The balance focus was 5.082 
inches aft of the 22.5° apex. Constructed of mahogany with an aluminum core, the 
drogue model did not decline from the horizontal when mounted. Sketches of the 


sting balance models are presented in Appendix D. 


25 


vntnn | eta EMBoe > 


eet tll ble 





Figure 3.12 Left Helix Sting Model 
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E. STATIC CONDITIONS SIMULATION MODEL 

An existing simulation model by Clifton was used to determine the Reynolds 
number range for wind tunnel testing [Ref.2]. The model described the static, no 
wind equilibrium conditions of a cable towed by an airplane in a circular orbit for 
positive angles of attack to 90°. Additionally, modifications were conducted by 
Clifton for inclusion of experimentally derived force coefficients for improved static 


predictions. 
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IV. EXPERIMENTAL PROCEDURES 


A. GENERAL 


1. Model Preparation 

Wall balance models required approximately two hours of preparation 
before nr: The steel cylindrical rods and support stanchions were sanded with 
fine grit finishing paper and soaked in an ethanol bath. After removal and drying 
they were coated in spray silicon and handled with latex gloves to prevent surface 
oxidation. The rods or wires were installed and set with two screws per piece per 
stanchion. ‘The set screw and end cap holes were filled with clay to reduce 
aerodynamic interference. The aluminum plate with model was then mounted on the 
balance. 

Sting balance models required less preparation. After surface cleaning 
with a soft cloth, the models were carefully slid over the Task balance and screws set 
into the locking recesses. Tape was used over the mounting holes on the drogue 
model and clay fill on the other sting models. 

2. Signal Conditioner / Amplifier Preparation 

The signal-conditioner span rheostat controlled bridge excitation voltage. 

For the wall balance channels, excitation was set at 10.00 Voc -0.0SmV (MC-MIO- 


16L-9 I/O board saturation was at 10 V,,). For the sting balance, channels N1 and 
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N2 excitation were set to 5.15 and 5.75 Vp. as required from the local calibration 
(Appendix B). The remaining four channels set 5.00 V,,. with all error +0.02mV. 

The Pacific’ amplifiers required input and output to be zeroed after 
model changes, long test periods and long off periods. With the amplifier input 
shorted and a gain of 1, the output set screw was adjusted to zero to +100 pV 
tolerance. After increasing the gain to 1000, the input screw was zeroed to +500 pV 
(Ref. 16]. Finally, the shorting plugs were pulled and the incoming signals set to 
zero, +0.01 mV with the signal control on the signal-conditioner panel. Acquisition 


system preparation was complete. 


B. WALL BALANCE EXPERIMENTS 


1. Test Matrix 

Simulation software by Clifton was run to determine the Reynolds number 
of the 20,000-foot-long 3x7 wire for test matrix construction [Ref. 2]. Results are 
presented in Figure 4.1. The simulation model used a density look-up table from 
flight test data, and constant viscosity. Hand calculations of Reynolds number using 
altitude varying temperature, density, and viscosity produced a low bound value of 
3500, and an upper bound Reynolds number of 17,500 [Ref. 3:p. 6-10]. After 
accounting for tunnel test conditions, test velocities were selected and are presented 


in Table 4.1. 
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Figure 4.1 Static Simulation - Reynolds Number vs Altitude 


TABLE 4.1 WALL BALANCE TEST VELOCITIES / REYNOLDS NUMBERS 


Reynolds Number 
3x7 Wire (d=0.158" Rods (d=0.3125" 


7 





The lowest dynamic pressure of 10 Ib,/ft? was selected due to balance resolution 


limitations. Velocities representative of Q=50 Ib,/ft? and Q=60 Ib,/ft? were not run 


30 


for the 0.3125-inch-diameter steel cylindrical rods to maintain Reynolds number 
ranges, and allow correlation of 3x7 wire data at similar dynamic pressures. All wall 
balance tunnel runs produced subcritical Reynolds numbers [Ref. 3:p. 106}. 

Wall balance tunnel experiments were started with the model at 105° to 
the flow (Figure 3.8) and advanced through 90° (Figure 3.9) to -15° in 5° 
decrements; resulting in 25 angle-of-attack settings. Four additional check points 
were taken per run. A total of 16 tunnel runs were completed in 11.4 tunnel 
operating hours. Breakdown of runs are as follows: a) clean stanchions alone - 6 
runs (Q=10-60 psf), b) 0.3125-inch-diameter rods - 4 runs (Q=10-40 psf), c) 3x7 
wires - 6 runs (Q= 10-60 psf). 

2. Balance Tare Values 

Wall balance channels exhibited a linear drift with respect to time 
requiring tare readings immediately prior and after tunnel operation. The tare 
difference was prorated assuming equal time between angle of attack positions by a 


QuickBasic program contained in Appendix C, and stored for data reduction. 


3. Tunnel Operation 
Wind tunnel operating procedures are summarized below: 
a) After setting desired velocity (cm H,0) in the tunnel manometer, 
the tunnel was started and brought up to and stabilized at desired dynamic pressure. 
b) When stable, 3 data points were taken. 


c) Model angle of attack was advanced to the next position. 
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d) Tunnel velocity was adjusted to maintain test Q. 
e) Steps (b) through (d) were repeated for all test points. 
f) The tunnel was shut down and elapsed time, temperature logged. 


A detailed description of tunnel operating procedures is presented in Reference 10. 


C. STING BALANCE EXPERIMENTS 


1. Test Matrix 

To maintain Reynolds number integrity with the wall balance experiments, 
the 15-times scale sting balance models would be required to be tested at velocities 
well below balance resolution. A mid-tunnel range dynamic pressure of 40 lb,/ft’ was 
selected to utilize the full calibration range, minimize turbulence, and achieve force 
measurements greater than 2% of balance resolution. Experimentation was also 
conducted at Q=20 lb,/ft? to determine velocity dependence. At Q=20 and 40 lb,/ft? 
the model Reynolds numbers were 1.63x10° and 2.31x10° respectively, still subcritical 
for a circular cylinder [Ref. 3]. 

Baseline model mounting position was with the two mounting screws in 
front, facing the front model access window. Test runs were repeated with mounting 
holes rotated to top, back, and bottom reference positions for each model to ensure 
helix represencnien of an infinite wire, and reveal asymmetries from nose and body 
construction. Additional test runs were conducted with a 0.5-inch-wide strip of 80 


and 220 grit sandpaper attached to the base of the nose designed to disrupt nose 
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generated asymmetric flow separation and vortex shedding at the higher angles of 
attack [Ref. 15]. 

Sting balance models were tested from 0° to 20° in 4° increments and 
from 20° to 50° in 5° increments. The angle-of-attack positions were repeated in the 
negative direction. Six check points were taken for a total of 29 angle-of-attack 
positions per experimental run. A total of 24 data runs were conducted in 11.5 
tunnel operating hours. A breakdown of runs, model and test conditions are 


presented in Table 4.2. 


TABLE 4.2 STING BALANCE MODEL TEST CONDITIONS 


| 2 (b/t) sition | _ Sandpaper 


None,80,220 
None,80,220 
None 
None 
None 
None 


Right Helix Front None, 220 


Top None, 220 
Back | None 
Bottom | None 
Front None 


Left Helix Front None, 220 
Top None, 220 
Back None 
Bottom None 
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2. Balance Tares 
The sting balance acquisition program recorded wind-off force and 
moment tare values to be subtracted from wind-on data points. The sting balance 


exhibited no drift. 


3. Tunnel Operation 
Wind tunnel operating procedures for sting balance experiments were 


identical to wall balance operating procedures. 


4. Flow Visualization 
String tufts flow visualization was conducted with the right helix model, 
front mounted, at a Q=20 lb,/ft?. The model, shown in Figure 4.2, had 36 1.75 inch 
black tufts taped in three bands positioned 4, 10, and 15.5 inches aft of the nose. 
Four tufts were attached on each Vs wrap of the helix with end tufts positioned close 
to the groove edges. One run of 0.9 hours was conducted with photographs taken 


at the same angle-of-attack positions as the data runs. 


D. DATA REDUCTION 


1. General 
Raw experimental data was transferred from the IBM PS/2° via 1.4 
megabyte, 3.5 inch disk to a 486-33 MHz computer for data reduction. Any 
corrections to inputed angles-of-attack values were completed using Microsoft DOS 
5.0 Editor prior to importing into QuattroPro 3.0 spreadsheet software. All data 


manipulation and graphing were conducted in QuattroPro 3.0. 
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Figure 4.2 Right Helix Model - Flow Visualization (Wind Off) 


2. Wall Balance Data 
Ten spreadsheets, one for each 3x7 wire and cylindrical rod velocity, were 

constructed for the wall balance experimentation. Each spreadsheet imported four 
files of 4 columns (axial force, axial moment, normal force, normal moment) by 75 
rows (AOA data). The four files required per spreadsheet were: a) the clean 
stanchion test run, b) its calculated tare variation, c) the wire or bar test run, and d) 
its representative tare calculations. Resultant normal and axial forces and moments 
were calculated using 

NET FORCES | _ [ (TEST FORCES/MOMENTS - MODEL TARE) (4,1) 

ET MOMENTS} - (CLEAN STANCHIONS - STANCHION TARE) | 
Computations were conducted in wire/cylinder vice balance coordinate system and 


defined as follows: 
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perpendicular to wire/cylinder 
parallel to wire/cylinder 


Normal 
Axial 


Nondimensional force coefficients were calculated for each angle of attack in the 


normal direction using 


eames (4.2) 
Q*L*d 





normal force coefficient 

normal force ([b,) 

test dynamic pressure (lb,/ft’) 
wire/cylinder cumulative length (ft) 
wire/cylinder diameter (ft) 


Where: C 


rl oe, Olea 
Mo oh wt 


Additionally, a normal force coefficient referenced to the normal velocity component 


was calculated using 


N 


ae (4.3) 
0.5*p *(Vn)? *L*d 


Cy, Vcomp ~ 


normal force coefficient 
normal force (Ib,) 


Where: C 


N, Vcomp 


Vn_ = velocity component in the normal direction (ft/s) 
p = atmospheric density (lb, s?/ft*) 

L = wire/cylinder cumulative length (ft) 

d = wire/cylinder diameter (ft) 


Axial force coefficients, C,, were calculated using (4.2) substituting axial force, A, for 


normal force, N. 


3. Sting Balance Data 
Four primary spreadsheets were constructed for the four sting balance 
models. Each test run file consisted of 6 columns (3 forces, 3 moments) and 84 rows 
(AOA data). The clean cylinder, right, and left helix model spreadsheets contained 


the four different mounting positions and different dynamic pressure data runs as 
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applicable. One spreadsheet contained all drogue data runs, and three separate 
spreadsheets were opened to reduce sandpaper modified model test data. 

For the clean cylinder and helix models, normal, axial, and side force 
coefficients (C,, C,, Cy) were calculated using the respective form of (4.2) with 
normal, axial, and side forces (N, A, S), L = cylinder reference length, and d = 
cylinder diameter. Additionally, nondimensional pitching and yawing moment 
coefficients were calculated using 


PITCHING/YAWING MOMENT (ftlb,) 
Q«L?«d 


Cmyy = (4.4) 


pitching, yawing moment coefficient 
dynamic pressure (lb,/ft’) 

length of cylinder section (ft) 
diameter of cylinder (ft) 


OW ee Qe 
i ot oi ith 


For the drogue, force and moment coefficients were calculated in (4.5) 
and (4.6) referenced to the base area and drogue length. Graphs of coefficients 


versus angle of attack were plotted. 


N, A 


S ee (4.5) 
“A Qx(n/4) 4d? 
PITCHING MOMENT (ft lb 
Cm, = eles SA hes tA bel (4.6) 
Q*(n/4) +d? *L 
Where: Cy, = normal, axial force coefficient 
Cm, = normal moment coefficient 
N,A = normal, axial force 
Q = dynamic pressure (lb,/ft’) 
d = drogue base diameter (ft) 
L = drogue length (ft) 


3) 


E. STATIC CONDITIONS SIMULATION MODEL 

An existing simulation model by Clifton [Ref.2] was modified with the 
experimentally derived normal, axial, and side force coefficients of the 3x7 wire. 
Drag and lift coefficients for the drogue and moment modifications were also 
included. The FORTRAN program, TAC17, was renamed TAC17A and is contained 
in Appendix E. Eight test runs using Microsoft FORTRAN 5.1 and TAC17 and 
TACI7A programs were conducted on a 486-33 MHz computer with conditions 
described in Tables 4.3 and 4.4. 

TABLE 4.3 STATIC SIMULATION TEST CONDITIONS 


— Altitude | Bank Name of 


| a | taci7 { 1s6_{ 18325 | 34° | Operational | 
|B | taci7 | is6_| 19325 | 34° | Level | 
18,325 

bp | Taci7A| 156 | 19325 | 34 








TABLE 4.4 SIMULATION TEST RUNS 


YES (3*C,) 
—s | B_ | No __| Previous Baseline __ 
6 | D_ | NO __| Comparison 
D 


Compare (Cyon) _| 
8 D YES (3*Cy) | Observed Trail Angle 
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Outputs from these programs included radial, theta, and z position, tension, angle of 
attack, true velocity, and Reynolds number for each grid point of the wire. The 64 
Output files were imported into one QuattroPro 3.0 spreadsheet for plotting and 


analysis. 
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V. RESULTS AND DISCUSSION 


A. HOERNER MODEL VERIFICATION 

Hoerner presented an empirical relationship for lift and drag coefficients of 
inclined cylinders at subcritical Reynolds numbers. His relationship was based on the 
premise that the "fluid-dynamic pressure forces of such bodies only correspond to the 
velocity component in the direction normal to their axis." [Ref. 6] From that 


relationship, for a circular cylinder, 


= a 


C, = Cp pasic *Sin’ a *+COSO (2.4) 





Where: C, = drag force coefficient 
C, = lift force coefficient 
Cy, = normal force coefficient (equation (4.2)) 
a = angle of attack 


Hoerner’s relationships are presented in Figure 2.3. Instead of defining an axial 
force coefficient, C,, Hoerner added a wC, = 0.02 skin friction term. By 
transforming a body axis system of C, and C, into C, and C,, Hoerner’s relationships 


can be written as equation (5.1) with C,=0. 
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a : {e exc (5.1) 
Gr c 1 f 


Using wall balance data for the 0.3125-inch-diameter clean cylinder and the 
relationships presented by Hoerner in (2.4) and (5.1), lift and drag coefficients were 
calculated. Plotted against Hoerner’s empirical findings, the Cp) and C, of the clean 
cylinders are presented in Figure 5.1. C, of the cylinders showed agreement with 
Hoerner except in the region of 90° angle of attack. The difference was attributed 
to Hoerner using Cp pasic = 1.1 while the Cy 99°, which is equal to Cp p,sico WaS calculated 
from (4.2) as 1.193. The C, of the experimental cylinders showed close agreement 
with Hoerner’s closed form relationships. 

Similar analysis was conducted using 3x7 wire wall balance data and is 
presented in Figure 5.2. Cy 99° of the wire, calculated from (4.2), was 1.09. Wire C, 
and C, showed fairly good agreement with Hoerner’s closed form-equations with the 
maximum deviations found in the middle range angles of attack. With verification 
of experimental data on the circular rods in accord with historical empirical 


relationships, definition of force coefficients in all three axes commenced. 


B. NORMAL FORCE COEFFICIENT 


1. Wall Balance Data 
Normal force coefficients for the wall balance data at each velocity were 
calculated using the appropriate form of (4.2). Cy 9. values for the clean cylinders 


and 3x7 wires plotted versus test dynamic pressure are presented in Figure 5.3. A 
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Figure §.2 C,, C,, Hoerner vs AOA 
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WALL, BALANCE MODEIS 
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Figure 5.3 Wall Balance C, vs Q 





slight variation of C, with respect to Q was observed. The highest values of C,, at 
Q=30 Ib,/ft’, were selected for further calculation with the following justifications: 
1) the force coefficient of a circular cylinder of finite length is less than that for a 
cylinder of infinite length [Ref.17], and 2) higher C, values produce higher wire 
tension (worse case condition) in the Clifton simulation model [Ref. 2]. 

Normal force coefficients calculated with respect to their normal velocity 
components using equation (4.3) are presented in Figures 5.4 and 5.5. The figures 
confirmed the dependence of C, on the velocity component as presented by Hoerner 
[Ref. 6]. Additionally, C, values were constant for angles of attack above 50° for the 


clean cylinder data and above 60° for the 3x7 wire data. Breakdown of correlation 
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Figure 5.4 Clean Cylinders C vs AOA 
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Figure 5.5 3x7 Wires Cy yeomp VS AOA 


Aa, 


increased as AOA decreased due to wake influences from the model stanchions. 
Wall balance data was considered unreliable below 20° AOA. 
Analytical fit of C, data was accomplished with emphasis on the 45° to 


105° AOA range. The clean cylinder normal coefficient data was expressed by 
iG (ce—C,.). ‘Sisco, = 1,193 (5.2) 
The curve fit is presented in Figure 5.6. The 3x7 wire data fit a sin’a relation at the 


higher AOAs but deviated low in the midrange angles of attack. After much study, 


the 3x7 wire data was fit by the following relationship: 


Ca) = Cy oo *Sin” o - Kxsin?(2a) (5.3) 
Where: Cyoo = 1.09 
K = 0.08 


Presented in Figure 5.7, correlation between data and equation (5.3) extended to 


AOAs from 105° down well below 45°. 


2. Sting Balance Data 
Normal force and moment coefficients for all runs of the clean cylinder, 
right, and left helix models were calculated using (4.2) and (4.4). C, showed no 
variation between test runs with different mounting positions for all three models as 
seen in Figures 5.8, 5.9, and 5.10. The C, versus AOA curves for the right and left 
helix models were identical, again indicating the dominance of normal force on long 
slender bodies [Ref. 18] regardless of groove direction. The use of sandpaper to trip 


and disrupt the flow had no influence on C,, values. Most importantly, C, showed 
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Eqn (5.2) 


E 
> 
! 
8 
E 
3) 
e 
° 


40 60 
ANGIE, OF ATTACK (deg) 


r) 


| 


= 


Fan (5.3) 





= 
= 


= 
os 


NORMAL COEFFICIENT (Cn-Vfree) 
= 


80 





40 60 
ANGLE OF ATTACK (deg) 


Figure 5.7 3x7 Wire C, vs AOA 
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CLEAN CYLINDER STING MODEL 
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Figure 5.9 Right Helix Sting Model C, vs AOA 
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Figure 5.10 Left Helix Sting Model C, vs AOA 


no dependence on dynamic pressure at the subcritical test Reynolds numbers as 
presented in Figures 5.11 and 5.12 for the clean cylinder and right helix models 
respectively. Q independence allowed for accurate integration of the sting and wall 
balance data. 

Integration of normal force data from the two dissimilar balance tests 
Started by using relationships (5.2) and (5.3) generated from wall balance data. As 
revealed in Figure 5.13, equation (5.2) only agreed with the clean cylinder sting 
model in the +20° AOA range. Sting C, values were significantly higher at higher 
AOAs. Analysis of sting Cm, and dC,/da data revealed the pitching moment center 
forward of the balance focal point suggesting ogive-nose influence on the normal 


force. C, for sting balance models was calculated using a reference length, L, equal 
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CLEAN CYLINDER MODEL — FRONT MOUNTED 
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Figure 5.12 Right Helix Sting Model, Q Influence on C,, 
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Fan (5.2) 
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Figure 5.13 Clean Cylinder Sting Model C, Analytical Fit Attempt 


to the cylinder length and did not account for nose reference area. A change in 
reference area would produce just a scaling constant application to (5.2). By noting 
the Cy 45- value of the sting balance was 25% larger than the correct value’ of Cy 4s- 


from equation (5.2), a scaling constant of 1.25 was applied to (5.2) to produce 


CAO) srg = 1.25*(1.193 *sin’a) (5.4) 


Equation (5.4) agreed with C, data for the clean cylinder as presented in Figure 5.14. 
The scaling constant was extended to the right helix tests, and applied to (5.3) to 


yield 


' From the wall balance (5.2): Cyas-) = 1.193*sin?(45°) = 0.5965 
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ClO) crvg = 1.25+*[1.09+sin?a - 0.08 +sin’(2a)] (5.5) 
Agreement of (5.5) with C,, data for the right helix model was confirmed in Figure 


5.15. 


bottom 
=e 


Eqn (5.4) 
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Figure 5.14 Clean Cylinder Sting Model C, Analytical Fit 
C. AXIAL FORCE COEFFICIENTS 





Axial force coefficients were calculated using the appropriate form of (4.2) and 
axial force data from the sting balance data. Shown in Figure 5.16, C, remained 


constant for low angles of attack. Over the test AOAs, C, can be summarized by 


0.017 |a|<45° 5.6) 


C (a) = 
AG | 0.0 Ja |>45° 


The magnitude of C, and range of AOAs was consistent with Hoerner [Ref. 6]. 
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Figure 5.16 Right Helix Sting Model C, vs AOA 





a2 


Axial force coefficients calculated from the wall balance experiments were 
considered erroneous. The tare procedure did not consider the aerodynamic 
interference at the stanchion/wire juncture. Experimental results by Roshko [Ref. 
19] revealed alterations in drag and vortex shedding frequency caused by a splitter 
plate in the wake of a cylinder. The wall balance models were assumed to have 
similar effects with stanchion/wire axial force actually less than forces from the 
stanchions alone. The consistent and expected value of C, from the sting balance 


experiments provided added justification for neglecting wall balance C, values. 


D. SIDE FORCE COEFFICIENTS 

Side force coefficients were calculated using the appropriate form of equation 
(4.2) for the clean cylinder, right, and left helix sting balance models. Presented in 
Figures 5.17, 5.18, and 5.19, C, showed agreement between +20° for the clean 
cylinder model, and +12° for the helix models with respect to screw mounting 
positions. Although the 220 grit sandpaper tended to "settle" the higher Cy values 
at higher AOAs for the clean cylinder, it displayed no effect on C, values for the 
helix models. This independence highlighted the strong aerodynamic influence of the 
helical grooves over possible asymmetric vortex shedding of the nose. As with 
normal force coefficients, Cy values were independent of test dynamic pressure. 
Additionally, resolution of yawing moment center, within 0.5 inches of balance focus, 


supports using cylinder length as reference for coefficient calculations. 
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Figure 5.18 Right Helix Sting Model C, vs AOA 
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Figure 5.19 Left Helix Sting Model C, vs AOA 





Determination of a side force relationship with respect to angle of attack for 
the helix model began with the baseline clean cylinder model. From Figure 5.17 for 
the clean cylinder, 


dC 
=] = -7.68x10~4 degs"' = -0.044 rad™! (5.7) 
da Joy 


The negative bias was caused by the model mounted -1.6° nose down (-S coordinate 
direction). The right helix model produced a larger negative slope, and the left helix 
a slight positive slope as presented in Figures 5.18 and 5.19 respectively. Subtracting 


(dC, /da)-y, produced a net slope for each helix model. 


a5 


= -0.0444 rad"! (5.8) 


ns 
Qu 
Sits: 
“ 
“eee” 
x 
| 


ey 0.0619 rad"! (5.9) 
Therefore, a side force, small in magnitude, was directly dependent on the helical 
groove direction. An analytical expression for C, was determined from Figure 5.18 
and (5.8). Selection of a sin(3q) relation resulted from an observed C, peak around 
+30° AOA and a tendency of Cy toward zero at +50° AOA. By multiplying (5.8) 


by cos(3q@) and integrating, Cy became 


— -0,0148 *sin(3«) |x | <60° (5.10) 
a) = ; 
rf 0.0 |x |>60° 


Equation (5.10) produced a good fit with experimental data in the +12° AOA range. 

The model side force was hypothesized to be attributed to circulation caused 
by the helical grooves. Figure 5.20 shows the right helix model at 20° AOA to the 
flow. Orientation of the grooves would direct the flow down and around the bottom 
of the model inducing circulation to create a negative side force, a result supported 
by data. Figure 5.21 shows the opposite effect with the left helix model at 20° AOA. 
The circulation theory was supported with flow-visualization tests. Figure 5.22, 
showing the backside of the right helix model at 20° AOA, clearly showed alignment 


of adjacent tufts to the grooves during tunnel operation. The alignment was 
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observed at lower angles of attack only, and tuft position was dominated by direct 


flow at higher angles of attack. 
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Figure 5.22 Right Helix Model in Flow Visualization, 20° AOA, Back Side 


5 


E. DROGUE MODEL DATA ANALYSIS 
Force and moment coefficients were calculated for the drogue sting balance 
model using equations (4.5) and (4.6). Again, data showed no variation with change 


in dynamic pressure. C, and C, were transformed into C, and C, using 


C os C 
L : cosa -sina ; N (5.11) 
C 5 sina cosa C - 
DROGUE DROGUE 


Shown in Figure 5.23, the drag coefficient showed little variation over the tested 





range of angles of attack. By noting C, = C, at a=0, 


Cc 


D 


= 0.53 (5.12) 


Resolution of pitching moment data produced an aerodynamic center referenced to 


the apex at 


8 = 0.7558 (5.13) 
L AC 


The aerodynamic center was invariant with angle of attack. The lift coefficient, 


shown in Figure 5.24, was expressed by 
C, procue ~ 0 Clee (rad) | |<45° (0.785 rad) (5.14) 


with C = 0.42 ata = 45°. The drogue model displayed stall conditions at angles 


L, max 


of attack greater than 60°. 


58 


032 SCALE DROGUE. MODEL 


Q=10 psf 
aes 


Q=20 psf 


-20 0 2() 40 
ANGLE OF ATTACK (degs) 


Figure 5.23 Drogue C, vs AOA 
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Figure 5.24 Drogue C, vs AOA 
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F, STATIC CONDITIONS SIMULATION MODEL 


1. Program Modifications 


The TAC17 static equilibrium program by Clifton [Ref. 2] was modified 


into TAC17A using experimentally derived coefficients for the wire and drogue. 


Specifically, C, and C, of the wire were rewritten as 


Cha) = 1.09*sin’a - 0.08*sin?(2a)  O0<a<90° (5.3) 
C(a) = 0.017 0" <a <45 (5.6) 
i 0.0 45° <a <90° 
and Cy was added with 
C (.) _ -0.0148 *sin(3 a) 0 <a <60 (5.10) 
J 0.0 60° <a <90° 


For the drogue, C,, was increased from 0.41 to 0.53. Previously used dC, /da = 2.0 


was replaced by the expression 


0.61 * a (rad) 0.0<a<0.785 rad (0° <a<45°) 
0.785<a<1.05 rad (45°<a<60°) 
1.05<a<1.57 rad (60°<a<90°) 


C,(a) = 0.42 (5.14) 
0.20 


Additionally, the aerodynamic center was moved from X/L = 0.7411 to X/L 


0.7558. 


60 


2. Simulation Tests 

Simulations test runs were conducted to verify and validate an 
approximate -15° (to the right looking aft) wire trail angle observed in flight test. 
The conditions of observation were 18,325 feet of altitude, 156 KEAS, and 0° bank 
angle, into the wind. Simulation was conducted at recorded altitude and airspeed, 
but 3° left bank angle due to convergence constraints in the numeric model. 
Running TAC17A with C, =0 produced a wire trail angle = -2°, calculated from the 
change in radial positions of the first 100 feet of wire from the airplane. With Cy 
activated, the trail angle increased to -7.6° or Aangle = -5.6° in the observed 
direction. TAC17A was run using 3*C\y, 1.e. 3 * equation (5.10), with a Aangle = - 
14.2°. A 3*C, simulation was considered accurate because of the following 
conditions: a) the flight-test-observed trail angle was quantified visually, not 
measured from instrumentation, b) the experimental model represented only 1.4 
inches and 0.75 revolutions of a 20,290 foot wire with 130,200 revolutions, and c) a 
finite-length circular cylinder produces a smaller force coefficient than one of infinite 
length [Ref. 17]. 

Comparison of simulation conditions at 3° bank angle was conducted with 
the TAC17A improvement showing some large differences with its TAC17 
predecessor. Of significant interest were the variation of wire radial position, angle 
of attack, and tension, presented in Figures 5.25, 5.26, and 5.27 respectively. The 
large influence of C, on radial position was dramatically evident, even with small 


amplitudes of C,. The changes in wire and drogue force coefficients doubled the 
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wire angle of attack increasing all force coefficients due to their sine relationships. 
The increase in coefficients increased the wire tension at the towplane approximately 
25%. The predicted static tension at the towplane was well within safety margins for 
avoiding wire failure [Ref. 2]. 

Simulation tests were also conducted in the operational conditions of 
18,325 feet altitude, 156 KEAS, and 34° bank angle. Comparison of simulation 
conditions again produced significant differences between the TAC17A and TAC17 
models. The experimental coefficients for the drogue, and C,, for the wire, produced 
the greatest changes as shown by wire radial in Figure 5.28 and angle of attack in 
Figure 5.29. Inclusion of Cy, and 3*C, showed little variation from the NO C, 
TACI17A simulation. The influence of C, has doubled the wire length displayed to 
AOAs greater than 60° increasing the C, force. The increase in forces again 
increased the wire tension at the towplane, from 800 lb, to 980 Ib, (22.5%). The 
major improvement from TAC17A’s changes came from the increased altitude range 
of the wire as shown in Figure 5.30. The increased z coordinate range increased the 
verticality from 55% to 73%. Verticality is defined as "the altitude of the towplane 
less the altitude of the drogue divided by the length of the wire" [Ref. 2]. The 
increased static verticality directly translates into increased antenna capability. 
Preliminary evaluation of experimental coefficients in the Clifton dynamic simulation 
model showed a decrease of deviation from flight test data from approximately 10% 


low to 5% high. 
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Figure 5.26 Cy, Influence on AOA 
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3x7 WIRE TENSION 
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Figure 5.28 Cy, Influence on Radial Position, Operational 
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Figure 5.30 Cy, Influence on Wire Z Position, Operational 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The experimentally derived force coefficients of a wire with a helical groove 
behaved similarly to historical relationships with respect to normal and axial forces. 


Specifically for the tested 3x7 wire 


Ca) = 1.09*sin?a - 0.08+sin’(2a)  0<a<90° 


0.017 0° <a <45° 
Cia) = 
0.0 45° <a <90° 


Discovered in this experimentation was a side force coefficient, small in magnitude 
and present at lower angles of attack. Resulting from groove induced flow 
circulation and dependent on groove direction, the side force coefficient for the 3x7 
wire was expressed as 
-0.0148 *sin(3«) 0° <a <60° 
Cy(«) i ° ° 
0.0 60° <a <90 
For the system drogue, the experimental drag coefficient and the aerodynamic 


center, referenced to the apex, were found to be invariant with angle of attack as 


C, = 0:53, A = 0.7558 


Finally, the lift coefficient can be expressed as 
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0.61 * a (rad) 0.0<sa<0.785 rad (0° <a<45°) 
C,(a) = 0.42 0.785<a<1.05 rad (45°<a<60°) 
0.20 1.05<a<1.57 rad (60°<a<90°) 


The wire force coefficients, coupled with the experimentally determined drogue 
coefficients, allowed for refinement of a previous static simulation model providing 
a closer relation between simulation prediction and flight test data. Additionally, 
inclusion of the experimental side force coefficient explained an observed 
phenomena, and revealed the impact the coefficient has over the "infinite" length of 
wire presented in this application. 

Further research with helically grooved infinite cylinders at subcritical Reynolds 
numbers 1s recommended. Specific recommendations include: 


@ the use of more sensitive balances with expected forces at least 50% of balance 
resolution vice 2%, 


@ models with more helical revolutions, 

® models with different helical angles, 

@ a detailed wake study, and 

® quantification of wire trail angles through instrumented flight test. 
Additionally, the computational challenge to represent the three dimensional forces 
on a wire should be undertaken to model generic helically grooved wires and reduce 


expensive flight test experimentation. 
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APPENDIX A: WALL BALANCE CALIBRATION 


An external, cylindrical, reflection-plane (wall) balance was used to measure 
forces and moments of the wire at the higher angles of attack. The flush mounted 
balance was designed with two pairs of strain gage bridges orthogonally mounted on 
flexure links separated vertically by 26.5 inches. Each bridge circuit had four active 
legs for automatic temperature compensation. The lower bridge A, and the upper 
bridge B, each output an axial and normal voltage (E,,, E,,, E,,, E,,). The balance 
coordinate system consisted of the axial direction parallel with tunnel walls and flow, 
and normal direction perpendicular to the tunnel walls while at 0° angle of attack. 

Calibration was conducted using the acquisition system with amplifier gain set 
at 1000, and MC-MIO-16L-9 board gain at 1, resulting in a 4.88 mV resolution. 
Calibration procedure consisted of loading a rig mounted on the aluminum plate in 
the normal and axial directions at two different heights. Figure A.1 and A.2 present 
the calibration rig loaded normally at a height equal to 7.75 inches above the tunnel 
floor. The pulley apparatus translated vertically to maintain level attitude of the 
cable. Prior to loading the rig, the amplifiers and signal conditioners were zeroed. 
The span control on the signal conditioner was set to 10 Vyc -0.05 Vp for each 
channel. Amplifiers were shorted and output and input zeroed at an amplifier gain 
of 1 and 1000, respectively. Finally, the shorting plugs were removed, the channels 


zeroed, and the acquisition program started. Weights, measured to 0.001 Ib, 
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accunacywere suspended from the rig. Measurements of E,,, E,,, E,,, E,, were 
recorded as weights were incremented and decremented. Care was taken to relieve 
the load on the cable when changing the weights to prevent a displayed balance 


hysteresis. 





Figure A.1 Wall Balance Calibration Rig 
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_ 


Figure A.2. Wall Balance Calibration Rig 





Four calibration runs were conducted, two in the normal and two in the axial 
direction, with the cable at 10.5625 and 7.75 inches of height referenced to the tunnel 
floor. The data from the four runs were imported into QuattroPro 3.0 for analysis. 
Plots of balance voltage versus calibration load are presented in Figures A.3, A.4, 
A.5, and A.6. The figures reveal the linearity expected from elastic loading, and the 
small interaction between channel bridges (cross-talk). Figures A.5 and A.6 also 
reveal the reduced sensitivity in the E,, channel due to replacing two legs of the 
bridge circuit on a questionable flexure link (i.e. possibly buckled) by constant 
reference resistance gages. Linear regression was conducted in QuattroPro to 
determine dAE/dload for each of the 16 conditions. The maximum standard 


deviation of the slopes was 1.5%. 
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Figure A.4 Calibration Voltage vs Normal Load, h=7.75 inches 
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Figure A.6 Calibration Voltage vs Axial Load, h=7.75 inches 
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The goal of the calibration was the calculation of the 4x4 calibration matrix, 


[K]. When [K] is post-multiplied by balance output voltages, forces and moments are 


produced. 
Fiz} {AXIAL FORCE 
E 
(Ks bat AXIAL MOMENT (A.1) 
E,,| || NORMAL FORCE 
E, NORMAL MOMENT 


Calculation of [K] came from inverting and post-multiplying the 4x4 dAE/dload 


matrix in 


/ / 
(, K,| |¢4E/dA dAE,JdA ddE,JdN doE,JdN 


Ky, Ky, K 

K,, Ky, Ko; Ky| |dAE,JdA dAE,JdA dAE,JdN dAE,JdN 
oe 

Ky, Ky Ky, Ky) |dAE JdA dAE.JdA dAE,/dN dAE_,JdN 


K, 
[Kar Ken Kes Kus dAE,/dA dAE,JdA dAdE,J/dN dAE,,JdN|  ‘A-2) 


1 1 0 0 

(a-b) (a-b)’ 0 0 

og a ae 
0 0 (a-b) (a-b)'} 


10.5625 inches 
7.75 inches 


Where: (a-b) = height at upper position 
(a-b) = height at lower position 


The 4x4 dAE/dload matrix was defined as follows: 


Column (1): slopes of voltages, axial load, h= 10.5625 inches 
Column (2): slopes of voltages, axial load, h=7.75 inches 
Column (3): slopes of voltages, normal load, h= 10.5625 inches 
Column (4): slopes of voltages, normal load, h=7.75 inches 
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The calibration matrix can be defined by four (2x2) partitioned elements as 


Kl K2 
K3 K4 


iK] = (A.3) 








where K1 and K4 determine the direct axial and normal forces and moments, and 
K2 and K3 contribute the cross channel interaction terms. After three full 


calibrations, [K] was determined as 


0.009292 -0.007686 -0.000053 -0.000209 
-0.033079 0.246046 0.007737 0.003644 (A.4) 

0.000063 -0.000417 0.009682 -0.004241 

| 0.002432 ~-0.006519 -0.033848 0.126897 


[K] = 
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APPENDIX B: STING BALANCE CALIBRATION 


A six degree of freedom, one inch diameter, Mark XIV, internal Task balance 
was used to measure forces and moments for models representing lower angles of 
attack of the wire, and for the drogue. The 4.9 inch long balance had eight 3500 
wheatstone bridges: two normal force (N1, N2), two side force (S1, S2), two axial 
force (Al, A2), and two rolling moment (11, 12). The two axial force, and two rolling 
moment bridges were each connected in parallel to produce 175) circuits, and 
therefore six electrical outputs: N1, N2, $1, $2, A, |. Maximum balance loads were 
400 lb, in the normal channels, 200 lb, in the side force channels, 100 lb; axially, and 
250 inlb, of rolling moment. On loan from NASA-Ames Research Center through 
the Navy-NASA Joint Institute of Aeronautics, the balance was calibrated to a 5 Vpc 
bridge excitation voltage by NASA-Ames personnel. Additional calibration was 
required to ensure proper force measurement, and moment focal point position prior 
to testing. 

Local calibration consisted of setting the six bridge excitations at 5 Vpc and 
evaluating loads and moments. The two normal force bridges, located +2.0 inches 
from the focal point, were loaded with an 11.2 lb, load. Evaluation of N1, N2, and 
total normal force output was conducted while loaded at the N1 bridge (+2 inches), 
the focal point, and the N2 bridge (-2 inches). For example, when the balance was 


loaded on the N1 position, an accurate calibration would indicate 11.2 lb, from the 
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N1 channel, 0 lb, from N2, and 11.2 lb, total. Similarly, if loaded at the focal point, 
N1 and N2 should both indicate 5.6 lb, with a total normal force, N = 11.2 lb, A 
reading high or low would indicate the direction of change required for the channel’s 
bridge excitation. By conducting several experiments at all three positions, the N1 
and N2 bridge excitation voltages were determined as 5.15 Vp and 5.75 Vp. 
respectively to yield proper loadings. 

The two side force bridges, 3.3 inches apart, were loaded in the same fashion 
with no change in excitation voltage required. Axial force was qualitatively assessed 


for magnitude and coordinate direction. 
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APPENDIX B, TAB I: STING BALANCE CALIBRATION CONSTANTS 


BALANCE CALIBRATION 
Start Date: 7-24-87 | Completion Date: 7-28-87 
NASA-Ames Calibration Laboratory 


Invoice No.: 440517 Kind: Force 
Pin No.: 3 Size: 1.0 
Make: Task 14C Rig No.:2 


| 400.00 | 400.00 | 350 | 0.1667 | 100. | 4625 | 


Capacity Max 9 x CAL | CAL 
(lb,) Load Gage | Shunt | ROG 


100.00 | 10000 | 175 | —~ | SoK | 4618 


200.00 200.00 | 350 | 0.1375 | 100.K | 4623 


21.0 fulb | 20.83 tub | 175 | ——_| 50K | 4623 










S1 3.1309E-02 | -3.8153E-08 | 3.2073E-02 | -8.9316E-09 | 0.263 | 0.132 


3.0366E-02 | -3.8607E-08 | 3.1167E-02 | -7.2517E-09 | 0.315 | 0.153 
| 1 |] 3.0885E-03 | 2.5672B-09 | 3.0908E-03 | -2.4769E-09 | 0.042 


Degree of Fit = 2 Accuracy = 15 Int. Degree of Fit = 2 
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INTERACTION COEFFICIENTS 


Negative Value 


N1/N2 5,8036E-03 -1,0257E-02 | 
N1/A 0.0000E-00 0.0000E-00 
N1/S1__| _ -4.1655E-03 4.5396E-03 


N1/S2 0.0000E-00 0.0000E-00 
N1/1 -5.8079E-02 4.4940E-02 


N2/N1 5.1 778E-02 
N2/A 4.4056E-03 
N2/S1 9.0385E-03 
N2/S2 0,0000E-00 
N2/I 6.1125E-02 
A/NI 2.1217E-03 
A/N2 -9.1524B-04 
A/S1 0.0000E-00 
A/S2 0.0000E-00 
A/l 9.7148E-02 
S1/N1 7.1275E-03 
S1/N2 | 0.0000E-00 | — 0.0000E-00 


2.6377E-01 6.2581E-02 
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Coefficient Negative Value 
| 1/N1 -3.5945E-04 | 
| 1/N2 0.0000E-00 
/A 0.0000E-00 
1/S1 0.0000E-00 
1/S2 0.0000E-00 


N1/N2? -7.9499E-07 
N1/A? 0.0000E-00 
N1/S?? 1.9670E-06 
N1/S2? 0.0000E-00 
N1/P 3,2320E-04 
N2/N?? -5,2897E-06 
N2/A? 1.0467E-05 
N2/S1? 4.8493E-07 
N2/S2? 0.0000E-00 
N2/P 1.1773E-03 
A/NY 4.2547E-06 
A/N2? -4.5946E-06 
A/S? 0.0000E-00 
A/S2 0.0000E-00 
A/P 7.5001E-04 
S1/NV2 1.2923E-05 
$1/N2? 0.0000E-00 
S1/A2 4.0345E-05 
$1/S2? 0.0000E-00 


S1/2 | -2.4592E-03 9.3969E-04 
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|_si/nz*_ | 0.0000E-00 | _8.6265-06 _ | 
| S2/si?_| 2.7825E-06 | -9.9830E-06 
‘8.0007E-04 


/N2 | -1.1512E- 0.0000E-00 
0.0000E-00 | — 0.0000E-00 
| ysr 0.0000E-00 | — 0.0000E-00 


5.1560E-08 0.0000E-00 


1/N1? 0.0000E-00 -1.5497E-07 
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A. 


APPENDIX C: ACQUISITION PROGRAMS 


WALL BALANCE PROGRAM 
MULTI. BAS 


This program was written and compiled using LabWindows and 
QuickBasic 4.5. (used "bc /o multi" to compile) It’s purpose 
is to read and convert voltages from four channels connected to 
the strain gauges on the tunnel wall balance. The voltages are 
converted to normal and axial forces and moments with respect to 
the balance. It was written and modified by LT Tom Stuart and 
LT Dean C. Schmidt, 24 July 92. 


Variables explained 


eaa = Strain gauge voltage at point A in Axial direction. 
eba = Strain gauge voltage at point B in Axial direction. 
ean = Strain gauge voltage at point A in Normal direction. 
ebn = Strain gauge voltage at point B in Normal direction. 


AX = Axial force 

Max = Axial moment 
NORM = Normal force 
Mnorm = Normal moment 


PKK KK KEK K KK KKK KKK KK KK KAKA KEKE KH KK KEK KK HH WKH HK KKH HHH HK HH KKH KKH HHH HEH KKK 


REM SINCLUDE: ‘'C:\LW\INCLUDE\LWSYSTEM. INC’ 
REM SINCLUDE: 'C:\LW\INCLUDE\GPIB. INC’ 
REM SINCLUDE: ‘C:\LW\INCLUDE\FORMATIO. INC’ 
REM SINCLUDE: ‘C:\LW\INCLUDE\GRAPHICS. INC’ 
REM SINCLUDE: 'C:\LW\INCLUDE\ANALYSIS. INC’ 
REM SINCLUDE: 'C:\LW\INCLUDE\DATAACQ. INC’ 
REM SINCLUDE: 'C:\LW\INCLUDE\RS232. INC’ 
DIM K#(4,4) 


DIM ean.array#(1000) ,eaa.array#(1000) ,ebn. array#(1000) ,eba. array#(1000) 
COMMON SHARED ean.array#() ,eaa.array#(),ebn.array#() ,eba.array#() 


DECLARE SUB volt (ean#, eaa#, ebn#, eba#) 


SCREEN 9, 0 
SOLOR 15, 1 


PRK KKA KKK KKK KK KKK KKK KKK KAKA HK KH HHH KKK HHH HHH HEH HE EHH HH HEH WHEE EEE HE Dk vk Ok vk Ok Dk ok 
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‘’ CALIBRATION MATRIX INPUT (See thesis for explaination) 


DATA 0.009292, -0.007686, -0.000053, -0.000209 
DATA -0.033079, 0.246046, 0.007737, 0.003644 
DATA 0.000063, -0.000417, 0.009682, -0.004241 
DATA 0.002432, -0.006519, -0.033848, 0.126897 


FOR L& = 1 TO 4: FOR M% = 1 TOG 
READ K#(L%,M%) : NEXT M% 
NEXT L% 


Paka ok de ak ok oe ak ok oe oe ak ok ak ok oe oe ake ak ak ok He ak He kak ak ke ek ak ak ak ak ak kok ok ak ak sk ok ak ak ak ak ak ak ok ak ok ak ak ak ok ok ak ak ok ae ak ok oe ok ok ak ok de ak ok ae ok 


LOCATE 10, 20: INPUT "Type the name of the voltage file"; VOLS 
VOLS = “C;\LW\ INSTR” 4 VOLS =e een 
OPEN VOLS FOR APPEND AS #1 


LOCATE 10, 20: INPUT "Type the name of the FORCE / MOMENT file"; FORS 
FORS = "C:\LW\INSTR\" + FORS + ".PRN" 
OPEN FORS FOR APPEND AS #2 


CLS: LOCATE 10, 20: PRINT "Make sure ‘Caps Lock’ is on now." 
SLEEP 2 


500 

CLS: LOCATE 10, 20: INPUT "Input the Test AOA"; alpha# 
alpha# = 90 - alpha# 

CLS: LOCATE 5, 20: INPUT “Continue? (Y/N)"; AS 


IF AS = "Y" THEN CALL volt (eaa#, eba#, ean#, ebn#) 
IF AS < "Y" THEN GOTO 5000 


PRINTS 

PRINT " AOA EAA (mV) EBA (mV) EAN (mV) 
EBN (mV)" 

PRINT " Shik kokbk kk ink bk Shik nbkkbkokk ak ok kak ok ak kok obok 
Shr kkk koko 

PRINT USING ” HHHHHH .HHHHHH": alpha#; eaa#; eba#; ean#; ebn# 

PRINT #1, USING "###### .######,"; alpha#; eaa#; eba#; ean#; ebn# 


‘ FORCES AND MOMENTS CALCULATIONS (See thesis for explaination) 


AX# = K#(1,1)*eaa# + K#(1,2)*eba# + K#(1,3)*ean# + K#(1,4)*ebn# 
Max# = K#(2,1)*eaa# + K#(2,2)*eba# + K#(2,3)*ean# + K#(2,4)*ebn# 
NORM# = K#(3,1)*eaa# + K#(3,2)*eba# + K#(3,3)*ean# + K#(3,4)*ebn# 
Mnorm# = K#(4,1)*eaa# + K#(4,2)*eba# + K#(4,3)*ean# + K#(4,4)*ebn# 


PROUN Ta = 


PRINT " AOA AXIAL (lb) MOMax (in-lb) NORMAL (lb) 
MOMnorm(in-lb)" 
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~ PRINT ” ok kk Skok obob bok ok okotot Shinn ook kok ok ok ok ok ook ok 
Sin bibsbot bk"! 


PRINT USING " HHHHH HHHHHH": alpha#; AX#; Max#; NORM#; Mnorm# 
PRINT #2, USING "##### .######,"; alpha#; AX#; Max#; NORM#; Mnorm# 


LOCATE 23, 15: INPUT "Do you want another reading? (Y/N)"; ANSS 
IF ANSS = "Y" THEN GOTO 500 

5000 CLOSE #1 

CLOSE #2 

END 


PSR HH KE HK HK He HH HH WH HHH HHH WH WH HHH WM WH HW HWW HHH BH WMH WH WH HHH HH HH WH HHH He WH Hee HE 


SUB volt (eaa#, eba#, ean#, ebn#) 
I lo nrin I ininiink ini ininbbiiinioiininivninbininioininigiboininikisbiirbiakaiotobsb sbi aibvaiink ak ok k ak ak aki oaks ok bak akok 


¢ 


' §$/R to read Channel 0,2,4,6 on MIO-16L-9 for Analog Voltage 


¢ 


Pk ak ok ok ok ok ok ok o& ve ok a ok Oe ak ok a ok oe Kk Hk ok ak ok Hk ke Kk ok ok ok Mk ke ek Mek ke Hk He ok Mek kk ak ok oe ok ae He ok ok He ok Fete ok de ok 


' Setting Board code for MIO-16L-9 
board. code%=0 


Bake ve ee oe ake eae OE ae HEME HE Me EE He WEE EH WWM BE WEME MW WH HH WM WM WM WE WME HE WEE HE HEM BW WE ME EEE WE MW ME sk ke ke Wek dhe gk 9k she 9k dk 


errl.num% = Init.DA.Brds(l, board.code%) 
err2.num% = AI.Setup(l, 0, 1) 
err3.num% = AI.Setup(l, 2, 1) 
err4.num% = AI.Setup(l, 4, 1) 
err5.num%’ = AI.Setup(l, 6, 1) 


' Configure and set clock to 1 MHZ 


err6.num% = CTR.Clock (1, 1, 1, 1) 
err/.num% = CTR.Config (1, 1, 0, 0, 0, 0) 


LWtotal! = 0 
FOR i = 1 TO 1000 
err8.num% = CTR.EvCount (1, 1, 1, O) 
‘ CHAN O = eaa 
err9 .num% = AI.Read(1l, 0, 1, value0$%) 


erl0O.num% = AI.Scale(1, 1, value0%, eaa.array#(is) ) 


' CHAN 2 = eba 
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erll.nums = AI 
erl2.num% = AI 


' CHAN 4 = ean 


erl3.num% = AIL. 


erl4.num% = AI 


' CHAN 6 = ebn 


erl5.num% = AI. 
.Scale(1l, 1, value6%, ebn.array#(i%) ) 


erl6.num% = AI 


erl7.num% = CTR.EvRead 


LWtotal! = LWto 


NEXT 1% 


.Read(1, 2, 1, value2%) 
.Scale(l, 1, value2%, eba.array#(it%) ) 


Read(1l, 4, 1, value4%) 


.Scale(1, 1, value4%, ean.array#(i%)) 


Read(l, 6, 1, value6%) 


(1, 1, overflo%, tcount$) 


tal! + tcount$ 


CLS: LOCATE 5,15:PRINT "Total Time is " LWtotal!*1lE-6" seconds." 


CALL Mean (eaa. 
CALL Mean (eba. 
CALL Mean (ean. 
CALL Mean (ebn. 


array#(), 1000, eaam#) 
array#(), 1000, ebam#) 
array#(), 1000, eanm#) 


array#(), 1000, ebnm#) 
Nokon iiinbiiniiniinkiiniiniininiinbihiiniinii i iii Ee aE eae ab ak ab ak ook 


’ This multiplication (*1000) will make the voltages in mV 
' This factor is required to work with the "K" cal matrix 


eaa#=eaam#* 1000 
eba#=ebam#* 1000 
ean#=eanm#* 1000 
ebn#=ebnm#* 1000 


END SUB 
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B. STING BALANCE PROGRAM 


: STING. BAS 


‘ This program was written and compiled using LabWindows and 

‘ QuickBasic 4.5. (used "bc /o multi" to compile) It's purpose 

‘ is to read and convert voltages from six sting balance channels 

‘ mounted in the Academic wind tunnel. The voltages are converted 

‘ using NASA-AMES balance calibration constants and equations written 
‘ by Yuan. The Labwindows program was written and modified by LT Tom 
‘ §tuart with assistance from LT Dean C. Schmidt. 

‘ Last date of modification: 10 Oct 92. 


‘ Variables explained 


f Nl = balance voltage at position 1 in the normal direction. 

; N2 = balance voltage at position 2 in the normal direction. 

: A = balance voltage in the axial direction. 

Sl = balance voltage at position 1 in the side force direction. 
: S2 = balance voltage at position 2 in the side force direction. 
: RM = balance voltage from rolling moment gauge. 


Perk kak wk ok HK KE ok ok ok a Ob Ok ok ak ak Ih Ok ak ok ok ak Mb ak Ok Hak ok ak ob ak hak ok Brak ak ok ak ob ae ak hak ak ok dhrok ok ok Oh ok ak ak ak OF ob ob ok ok ok ok Hob ok 


REM SINCLUDE: ‘C:\LW\INCLUDE\LWSYSTEM. INC’ 

REM SINCLUDE: 'C:\LW\INCLUDE\GPIB. INC’ 

REM SINCLUDE: 'C:\LW\INCLUDE\FORMATIO. INC’ 

REM SINCLUDE: ‘C:\LW\INCLUDE\GRAPHICS. INC’ 

REM SINCLUDE: ‘C:\LW\INCLUDE\ANALYSIS. INC’ 

REM SINCLUDE: ‘C:\LW\INCLUDE\DATAACQ. INC’ 

REM SINCLUDE: 'C:\LW\INCLUDE\RS232. INC’ 

D I M 


none. array#(1000) ,ntwo.array#(1000) , axial.array#(1000) ,sone.array#(1000) 
DIM stwo.array#(1000) ,rm.array#(1000) 

COMMON SHARED none.array#() ,ntwo.array#(),axial.array#() ,sone.array#() 
COMMON SHARED stwo.array#() ,rm.array#() 


DECLARE SUB volt (none#,ntwo#, axial#, sone#, stwo#, rm#) 


SCREEN 9, 0 
COLOR 15, 1 


Dokki iin niniibiinbi ok iinbiinbiniibiinintinninbiinbiinbin iin bb bb bbb 
‘ CALIBRATION CONSTANTS (See thesis for list of constants) 


‘ The calibration constants for the direct force nonlinear equations 
‘ and the force interaction equations are listed under separate appendix. 
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P ske ak ak ak ok ak ak ok ok ok ok ak ak HK ak ae a ok EK a HK ae ME HE EH ake ME HAH HEN EHH HAHAH HM WEE HWW HH HWW HH HW HH HHH HK WHA KK HK KH 


LOCATE 10, 20: INPUT "Type the name of the raw data file"; VOLS 
VOLS = "C:\LW\INSTR\" + VOLS + ".PRN” 
OPEN VOLS FOR APPEND AS #1 


LOCATE 10, 20: INPUT "Type the name of the FORCE / MOMENT file"; FORS 
FORS = "C:\LW\INSTR\" + FORS + ".PRN" 
OPEN FORS FOR APPEND AS #2 


CLS: LOCATE 10, 20: PRINT "Make sure ‘Caps Lock’ is on now." 
SEBER a2 


500 
CLS: LOCATE 10, 20: INPUT "Input the Test AOA"; alpha# 
alpha# = 90 - alpha# 


CLS: LOCATE 5, 20: INPUT "Continue? (Y/N)"; AS 
IF AS <> "N" THEN CALL volt (none#,ntwo#, axial#, sone#, stwo#, rm#) 
IF AS = "N" THEN GOTO 5000 


‘PRINT #1, USING "###### .#####4,"; alpha#; none#; ntwo#; axial#; sone#; 
Sstwo#; rm# 


1 feo ak ak ok ok ok ak ak ok ake ak ak ak ak ak ai ok ak ak ak ok ak ok ak ok ak ak ok ak ok ok ak ak ok ok 2b ak a 2k 2k ak ok ak ab ok aro ok ak ak ok ok ak ak ok ak ok ok ok ok ok ok ok 2k ok ak ak ak ab ok 
dkok ok ak ok 

DK EAA EAE EEK DE EEE IEE ETEK TE IEE IEE IE IIE IE aE aI aK a a a aaa ai aa ak a ok ak ak ak ak ak ok ak a ak ak ak ok ok ak ok ab ak ak ok ab ok 
deb okok 

' FORCES AND MOMENTS CALCULATIONS (Positive and Negative Equations) 

in niin niininininiiinininiiniinbinkinbiiinbiiniiniink bib ibibo bok akira ak abi 2k ak kok ak ak ak ak ab ok ak ak ok ak ok ak ok ab ok 
dkok bake 

VEX = 5! ‘Excitation voltage 

Nl = none#*5000! /VEX 

N2 = ntwo#*5000! /VEX 

A = axial#*5000! /VEX 

Sl = sone#*5000! /VEX 

S2 = stwo#*5000! /VEX 

R = rm#*416.6/7#/VEX 

600 IF none#>=0 THEN GOTO 1000 ELSE GOTO 2000 

620 IF ntwo#>=0 THEN GOTO 1100 ELSE GOTO 2100 

640 IF axial#>=0 THEN GOTO 1200 ELSE GOTO 2200 

660 IF sone#>=0 THEN GOTO 1300 ELSE GOTO 2300 

680 IF stwo#>=0 THEN GOTO 1400 ELSE GOTO 2400 

700 IF rm#>=0 THEN GOTO 1500 ELSE GOTO 2500 
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Mek POSITIVE EQUATIONS *** keke x 


¢ 


1000 EN1 = 0.050861*N1-5.4826E-09*(N1°2) 


GOTO 620 

1100 EN2 = 0.047211*N2-1.7015E-08*(N2°2) 
GOTO 640 

1200 EA = 0.014309*A-7.1962E-10*(A“2) 
GOTO 660 

1300 ES1 = 0.031309*S1-3.8253E-08*(S1°2) 
GOTO 680 

1400 ES2 = 0.030366*S2-3.8607E-08*(S2°2) 
GOTO 7/00 

1500 ER = 0.0030885*R+2 .5672E-09*(R“2) 
GOTO 2600 


’ 


'xkeeKKK NEGATIVE EQUATIONS ****xx** 

2000 EN1 = 0.051591*N1+1.7157E-08*(N1°2) 
GOTO 620 

2100 EN2 = 0.047763*N2+8.915299E-09%*(N2°2) 
GOTO 640 

2200 EA = 0.01429*A-1.3322E-09*(A"2) 

GOTO 660 

2300 ES1 = 0.032073*S1-8.931601E-09*(S1%2) 
GOTO 680 

2400 ES2 = 0.031167*S2-7.2517E-09*(S2°2) 
GOTO 700 

2500 ER = 0.0030908*R-2.4769E-09*(R*2) 


¢ 


'xkkkxe*eK* FORCE INTERACTION EQUATIONS sk kkk kk kok kk kk kk ek kt 

2600 

ERINT " " 

PRINT " eke kee FORCE INTERACTION CONVERGENCE **x%*%xx% 
PRINT " CYCLE AOA N1 N2 A Sl S2 

R" 

PRINT " # deg lb lb lb lb lb 
help " 

r R I N le 
Werke te kk ok kk ok ok kok ok kt ok kkk tok kk Kk tk kkk ke tk ke kkk ok kkk bok bok tk kkk ok ok tk ok ok ok kk ok ok ok ok ok 
di. i i a a 


; Iteration to check for convergence 

CYCLE = 0 

FOR I = 1 TO 10 

2800 IF none#>=0 THEN GOTO 3000 ELSE GOTO 4000 

2820 IF ntwo#>=0 THEN GOTO 3100 ELSE GOTO 4100 

2840 IF axial#>=0 THEN GOTO 3200 ELSE GOTO 4200 
2860 IF sone#>=0 THEN GOTO 3300 ELSE GOTO 4300 
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2880 IF stwo#>=0 THEN GOTO 3400 ELSE GOTO 4400 
2900 IF rm#>=0 THEN GOTO 3500 ELSE GOTO 4500 


i kikkkeK POSITIVE EQUATIONS **kkkkK 

3000 XN1 = EN1+0.0058036*N2+0.0041655*S1+0.058079*R-7.1926E-07*(N2° 2) 
XN1=XN1+4 .0352E-06*(S1°2) -0.0006786*(R*2) 

GOTO 2820 


3100 XN2 = EN2+0.046218%*N1-0.0028393*A-0.0081694%*S1+0 .0041463%*S2 


XN2=XN2+0 .077279*R-6.8577E-07*(N1°2)-1.7755E-05*(A*2)+2.1719E-06*(S1%2) 
XN2=XN241 . 8582E-06*(S2%2)-0.0019294*(R*2) 
GOTO 2840 


3200 XA = EA+8.6893E-044N1+6 .0359E-04%*S147.7722E-05*S2-0.11115*R 
XA=XA+4 .4537E-07*(N1°2)+4.7936E-06*(S1°2) -4.1033E-06*(S2°%2) 
XA=XA+2 .0967E-04*(R* 2) 

GOTO 2860 


3300 XS1 = ES1-6.3459E-04*N1-0.11148*R+5 .535E-06*(N1°2)+0.0024592*(R*2) 
GOTO 2880 


3400 XS2 = ES2-0.0024237*N14+0.0022455*A+0 .0066785*S1-0.26377*R 
XS2=XS24+1.7099E-06*(N1°2)+1.2072E-05*(A%2) -2.7825E-06*(S1°%2) 
XS2=XS2+0.0062217*(R° 2) 

GOTO 2900 


3 5 0 0 X R 
ER-1.9928E-04*N2-2.5893E-04%S241.1512E-07*(N2°2) -5.156E-08*(S2°2) 
GOTO 4600 


Nnkinkbkk NEGATIVE EQUATIONS *ixe kee 

4000 XN1 = EN14+0.010257*N2-0.0045396*S1-0.04494*R+7 .9499E-07*(N2°2) 
XN1=XN1-1.967E-06*(S1%*2) -0.0003232*(R*2) 

GOTO 2820 


4100 XN2 = EN2+0.051778*N1-0.0044056*A-9.038499E-03%*S1-0.061125*R 
XN2=XN2+5.2897E-06%*(N1°2)+1.0467E-05*(A%2) -4.8493E-07*(S1°2) 
XN2=XN2-0.0011773*(R* 2) 

GOTO 2840 


4200 XA = EA-0.0021217*N1+9.1524E-04*N2-0.097148*R-4.2547E-06*(N1°2) 
XA=XA+4 .5846E-06%*(N2°2)-7.5001E-04*(R*2) 
GOTO 2860 


4300 XS1 = ES1-0.0071275*N1-0.0089235*A-0.05268*R-1.2923E-05*(N1°2) 


XS1=XS1-4.0345E-05*(A%2) -9.3969E-04%*(R* 2) 
GOTO 2880 
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4400 XS2 = ES2-0.0037176*N1-0.0052619*N2+0 .0072915*A+0 .006856*S1 
XS2=XS2-0.062581*R-5.211E-07*(N1°2) -8.6265E-06*(N2°2)+3.7054E-05*(A*2) 
XS2=XS2+9 .983001E-06*(S1°2)+8.0007E-04*(R° 2) 
GOTO 2900 


4500 XR = ER+3.5945E-04*N1+1.5497E-07*(N2° 2) 


4600 ' Rename for next iteration 
Nl = XN1 

N2 = XN2 

A = XA 

Sl = XS1l 

S2 = XS2 

R = XR 


s 


: Counter for convergence iterations 
eyYCLE@=— CYCLE + 1 


PRINT USING “####.####"; CYCLE; alpha#; Nl; N2; A; Sl; S2; R 
NEXT I 
PRINT #1, USING "####.######"; CYCLE; alpha#; Nl; N2; A; Sl; S2; R 


LOCATE 21,15: INPUT "IS CONVERGENCE OK? (Y or N)";CONVS 
IF CONVS = "N" THEN GOTO 2600 


NORMAL = Nl + N2 

SIDE = Sl + S2 

AXIAL = A 

PITCH = (N1-N2)*0.1667 
YAW = (S1-$2)*0.1375 
ROLL = R/12.0 


TKK KK KKK KKK HHH HE HHH HE HEE EE EE IETF AEE III IE IE aE aE ab ab aE EEE ara ak a ab ak ok ak ak 
eek TARE CALCULATIONS ciiininkiiininnknkvkiiinkinknkibbbk kok ab kak ak ak ak ak abak ak ak abak 


LOCATE 23,15:INPUT "IS THIS A TARE READING? (Y or N)";TARS 
IF TARS<"Y" GOTO 4700 


TNORM = NORMAL 
io DE SIDE 
TAXIAL = AXIAL 
TPITCH = PITCH 
TYAW = YAW 
TROLL = ROLL 
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CLS 


PRINT 

PRINT " KkekKKKK TARE CALCULATIONS keke" 

PRINT " NORMAL SIDE AXTAL PITCH YAW 
ROLL" 

PRINT " lb lb lb ftlb ftlb 
ftlb 

P R I N T sh 


KKK KKK XK * kk Me ke HK Kk Hk ek ek Me ke ek ee HK KE ke kk ke kk kw ke ok kk ok ke ek ek ek ok ek ok kk ke ek KK 
KK Kk 

PRINT USING “#####. ######,": TNORM; TSIDE; TAXIAL; TPITCH; TYAW; TROLL 
PRINT #2, USING "#####.######,": alpha#; TNORM; TSIDE; TAXIAL; TPITCH; 
TYAW: TROLL 


GOTO 4800 
Np inininininniininni nnn nin rnin rrr rrr ae aided airs ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak oak ak ak ok ok abot ak ok 


HH RK HK WX KH KKK HH 


4700 ' 

kkk FORCE CALCULATIONS *x* xxx KKK 
NORMF = NORMAL- TNORM 

SIDEF = SIDE-TSIDE 

AXTIALF = AXIAL-TAXIAL 

PITCHF = PITCH-TPITCH 

YAWF = YAW-TYAW 

ROLLF = ROLL-TROLL 


Gus 

PRINT " " 

PRINT " Soke FORCE CALCULATIONS seeibek 

PRINT " AOA="; alpha# 

PRINT " " 

PRINT " NORMAL SIDE AXIAL PITCH YAW 
ROLL" 

PRINT " 1b Lb lb ftlb ftlb 
ftlb 

P R I N iL " 


FHKE HAR HELE HH HEHE HA HEHEHE HHH HEHE FEE ETE HHH HHH HHH WEE HHH FEE TE TE EE ETE TE EH HEHE TE TE TE TE TE ETE DEDEDE 
ok k bok" 

PRINT USING "#####. ######,"; NORMF; SIDEF; AXIALF; PITCHF; YAWF; ROLLF 
PRINT #2, USING "“#####.######,"; alpha#; NORMF; SIDEF; AXIALF; PITCHF; 
YAWF; ROLLF 


4800 LOCATE 23, 15: INPUT "Do you want another reading? (Y/N)"; ANSS 
IF ANSS <> "N" THEN GOTO 500 

5000 CLOSE #1 

CLOSE #2 

END 
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kkk ke kok kk KK kK Kk Hk Kk Hk HH HAHAHAHAHA AHHH HK HHH HHA HH HH HHH HAN HHH HHH HW HHH WH HK WH He 


SUB volt (none#,ntwo#, axial#,sone#,stwo#, rm#) 
Pen nr 


U 


‘ S/R to read Channel 0,1,2,3,4,5 on MIO-16L-9 for Analog Voltage 


é 


kok kako Kk kK Kk Kk HH HH HH HHH HAHAHAHAHA AHHH AHHH HHH KKK ae eK HK eK sek Hak ak vk ak HK ve HK ek HK ae vk Hk 


‘ Setting Board code for MIO-16L-9 
board. code%=0 


BKK K HHH HHH HHH HAHAH HH HAHAH NHHH HH HH AK HHH HH HH HHH HHH HH HH HHH HK HHH NH HWE KKK HK HK HH WK KKK 


errl.num% = Init.DA.Brds(l, board.codes) 
err2.num%’ = Al.Setup(l, 0, 1) 


err3.num% = Al.Setup(l, 1, 1) 
err4.num% = Al.Setup(l, 2, 1) 
err5.num% = AlI.Setup(l1, 3, 1) 
err6.num% = Al.Setup(l, 4, 1) 
err/7.num% = Al.Setup(l, 5, 1) 


’ Configure and set clock to 1 MHZ 


err8.num% = CTR.Clock (1, 1, 1, 1) 
err9.num% = CTR.Config (1, 1, 0, 0, O, QO) 


LWtotal! = 0 

FOR i% = 1 TO 1000 

erl10.num% = CTR.EvCount (1, 1, 1, 0) 
‘ CHAN O = none 


erll.num% = AI.Read(1, 0, 1, value0%) 
erl2.num% = AI.Scale(1l, 1, value0%, none.array#(i%)) 


‘ CHAN 1 = ntwo 
erl3.num% = AI.Read(l, 1, 1, valuel%) 
erl4.num% = Al.Scale(1, 1, valuel%, ntwo.array#(i%) ) 


‘ CHAN 2 = axial 
erl5.num% = AI.Read(l, 2, 1, value2$%) 
erl6.num% = AI.Scale(1, 1, value2%, axial.array#(i%)) 


‘ CHAN 3 = sone 


erl7.num% = AlI.Read(1l, 3, 1, value3%) 
erl8.num% = AI.Scale(1, 1, value3%, sone.array#(i) ) 
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‘ CHAN 4 = stwo 
erl9.num% = AI.Read(1, 4, 1, value4$%) 
er20.num% = AlI.Scale(1, 1, value4%, stwo.array#(i$) ) 


’ CHAN 5 = 


rolling moment 


er21.num% = AI.Read(1, 5, 1, value5$) 
er22.num% = AlI.Scale(1, 1, value5%, rm.array#(i$) ) 


er23.num% = CTR.EvRead (1, 1, overflo%, tcount$%) 


LWtotal! = 


NEXT 1% 


LWtotal! + tcount% 


CLS: LOCATE 2,15:PRINT "Total Time is " LWtotal!*1lE-6" seconds." 


CALL Mean 
CALL Mean 
CALL Mean 
CALL Mean 
CALL Mean 
CALL Mean 


(none.array#(), 1000, none#) 
(ntwo.array#(), 1000, ntwo#) 
(axial.array#(), 1000, axial#) 
(sone.array#(), 1000, sone#) 
(stwo.array#(), 1000, stwo#) 
(rm.array#(), 1000, rm#) 


PKK HKHKK HH HHK HHH HH HHH KKK HK WK KKH HH HHH WKH HHH HAH HHH HH HH HHH HK HK Hk Hk vk vk ak kk 9k sk vk sk vk ak yk 


END SUB 
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C. WALL BALANCE TARE CALCULATION PROGRAM 


; TARE CONVERSION CALCULATION 


' The purpose of this program is to calculate the time varing 


‘windoff tare readings given the four channel tare values immediately 


‘prior to tunnel operation, and immediately after shutdown. 
‘Assumptions include: linear drift, samples taken at equal time 
‘intervals. This program is designed to run in Quickbasic 4.5. 


Dare AC4. 4), BCS), CC35, 4) 
TARES = "C:\QPRO\XXX. PRN" 
OPEN TARES FOR OUTPUT AS #1 


P = 31 ‘Total number of data points taken 
T= .4 ‘Total run time 
delT = T / (P - 1)’Incremented time step per point 


DATA -0.722499, 16.772749, 0.475904, -9.318065 
DATA -1.436443, 39.308192, 0.403325, -5.661456 


FOR L = 1 TO 2 
FOR M = 1 TO 4 

READ A(L, M) 
NEXT M 

NEXT L 


FOR N = 1 TO 4 
B(N) = (A(2, N) - A(1, N)) / T ‘Calculates the slope 
NEXT N 


AOA = 105 
FOR X = 1 TO P - l 


FOR K = 1 TO 4 
C(X, K) = A(1l, K) + BC(K) * delT * X ‘THE Calculation 
NEXT K 
nOKeS = | TO 3 
PRINT #1, USING "“##### .#####4#,"; AOA; C(X, 1); CCX, 2): CCX, 


(XxX, 4) 
NEXT S ‘The values are printed three times to match the three 
‘samples taken at each AOA during the runs 
AOA = AOA - 5 ‘Decrement AOA 
NEXT X 
CLOSE #1 
END 


95 


Os 


S9YIIIYS [OPO aouLyeg [[PAA T° ans14y 





SIXY SONVIVE AWM | ae oe 
M3IA JNOUS NOILVLOY OL sin0a 
P 3 ivId SNOIGN3IWIG 
= asve waiv, “3!'* Nwid 
pe aan 


“2 ) 
Se = ie io as 


AGIM VNNSLNY 


- ~9SIO (S 
SOOM WWNaG PG ee 28 (L*E)— 9) 





ONIQVdS O02 IV SQ0Y% 1140 
SUBONINAD 00°21-(9) AALS .g\E°O ty) 
“SDOIANOD YSQANIIAD 


APPENDIX D: SKETCHES OF EXPERIMENTAL MODELS 


1BBLS 5EVO*SEI 


96 


SaYIIOAS [aPOW souvyeg sung sapuyAD 7q ainsiy 


MNV_ESG SNOG3» 


30 S\ =319NYV XIV3N 2 
QaovasgnlL 3ididi 34v 





BACOAD SaA0COYD AWI113nh | 
Hid30 : 
1S2°0 | ae 








L8°S2 
(ADN3ZDNVL 3AID0) O9712 
0012 7 


ale 


 RAVAWD —— 
BONVIVA SNS A 
NO1D3a aN0o0d» AvDIISH 


WOANINOAAW 2 TaLVW 2 
SAHINI NI SNOISNAWIO ‘| 


97 


snjeieddy aouvyeg sung p’q ainsiy 


JePpo|] aouBleg suns ansoiq ¢g amnsiy 








SIxy | SISVINSOL 
SAHNI Bav O1 .108 
SNOISNBWIO NOlLWLON FP Boor 
[ corccrwsese ares ae Me 
o-2 502-4 | | 
ay aaa 
ae Oe sO | 
600"! O'g2 
sazos] .. 
3onvive} °"° A 
O'2 ey ~666'0 | : oe 
eee 
is ' <oeo 
O06'b s 
oo! | 3A331S 


"WO TY 


ANWVSONVW 
S3BHON\ 3YV 


SNOISNAWIO 








Ove S‘O gO? | 


meas | vIideckidikhiit 


ONin13a9~” 


APPENDIX E: STATIC CONDITIONS SIMULATION MODEL 


PROGRAM TAC1/7A 

THIS PROGRAM IS AN ITERATIVE SOLUTION TO THE STEADY 

STATE TACAMO WIRE PROBLEM. 

THE COEFFICIENTS AND DIMENSIONS REFLECT THE NEW 3X7 RATTAIL 


WIRE. 
CHEK HAHAHAHAHAHA HAHAHAHAHA HHH HAHAHAHAHAHA BEBE EH DE TEEPE TE TEE DEDEDE TE TE DE TE EE DE TE KH IK 


Ca Cane) 


C DECLARE AND DIMENSION VARIABLES. 
e 
C SCALARS FIRST. 
G 
é Al, Bl, Cl, A2, B2, C2 ARE THE TENSION TIMES THE SLOPE AT 
€ THE HALF STEP POINTS USED IN THE ITERATIVE SOLUTION. 
REAL Al,B1,Cl 
REAL A2,B2,C2 
C A22, B22, C22 ARE PLACEHOLDERS USED DURING THE ITERATIONS. 
REAL A22,B22,C22 
C A2A2, B2B2, C2C2 ARE SUMMERS USED IN AVERAGING. 
REAL A2A2,B2B2,C2C2 
6 AC IS THE AERODYNAMIC CENTER OF THE DROGUE MEASURED FROM THE LE. 
REAL AC 
C ALFAD IS THE DROGUE ANGLE OF ATTACK. 
REAL ALFAD 
C ALTTP IS THE TOWPLANE ALTITUDE IN FEET. 
REAL ALTTP 
© ASFD IS THE AERODYNAMIC SIDEFORCE OF THE DROGUE DUE TO THE 
6 SIDESLIP ANGLE BETA. 
REAL ASFD 
C BETA IS THE SIDESLIP ANGLE OF THE DROGUE. 
REAL BETA 
C CDD IS THE CD FOR THE DROGUE. 
REAL CDD 
6 CA IS THE WIRE AXIAL AERODYNAMIC FORCE COEFFICIENT. 
REAL CA 
6 CN IS THE WIRE NORMAL AERODYNAMIC FORCE COEFFICIENT. 
REAL CN 
C CG IS THE CENTER OF GRAVITY OF THE DROGUE, MEASURED AFT OF THE 
C Ties 
REAL CG 
G CLALD IS THE LIFT COEFFICIENT CURVE SLOPE FOR THE DROGUE. 
REAL CLALD 
C CLALDM IS THE MAXIMUM LIFT CURVE SLOPE FOR THE DROGUE. 
REAL CLALDM 
C CMACD IS THE COEFFICIENT OF MOMENT AROUND THE AERODYNAMIC 
C CENTER FOR THE DROGUE. 
REAL CMACD 


9° 


CY IS THE WIRE SIDEFORCE COEFFICIENT. 

REAL CY 

D IS THE WIRE DIAMETER. 

REAL D 

DELTAS IS THE INCREMENT OF WIRE LENGTH AT THE N‘'TH GRIDPOINT. 
REAL DELTAS 

G IS THE ACCELERATION DUE TO GRAVITY. 


REAL G 
LEND IS THE LENGTH OF THE DROGUE. 

REAL LEND 

LD IS THE LIFT PRODUCED BY THE DROGUE. 
REAL LD 

MHU IS THE MASS OF THE WIRE PER UNIT LENGTH. 
REAL MHU 

PHI IS THE ANGLE OF BANK OF THE TOWPLANE. 
REAL PHI 

Q IS THE LOCAL DYNAMIC PRESSURE. 

REAL Q 

RADTP IS THE ORBIT RADIUS OF THE TOWPLANE. 
REAL RADTP 

RHO IS THE LOCAL AIR DENSITY. 

REAL RHO 


RNEW IS A PLACEHOLDER FOR R(1) WHILE AN UPDATE IS CALCULATED 
ON THE OUTER LOOP. 

REAL RNEW 

STUFF1,2,3,4,5,6,7 ARE DUMMY VARIABLES FOR INTERIM CALCULATIONS. 
REAL STUFF1, STUFF2, STUFF3 

REAL STUFF4, STUFF5, STUFF6 

REAL STUFF/ 

THEDOT IS THE ORBIT RATE IN RADIANS PER SECOND. 

REAL THEDOT 

RR IS A CONSTANT USED TO START THE RADIAL COORDINATE 
CALCULATIONS AT THE DROGUE. 

REAL RR 

SD IS THE MAXIMUM CROSS SECTIONAL AREA OF THE DROGUE. 
REAL SD 

THTH IS A CONSTANT USED TO START THE THETA COORDINATE 
CALCULATIONS AT THE DROGUE. 

REAL THTH 

VEAS IS THE TOWPLANE EQUIVALENT AIRSPEED. 

REAL VEAS 

VTRUE IS THE TOWPLANE TRUE AIRSPEED. 

REAL VTRUE 

WD IS THE WEIGHT OF THE DROGUE. 

REAL WD 

ZNEW IS A PLACEHOLDER FOR Z(1) WHILE AN UPDATE IS 
CALCULATED ON THE OUTER LOOP. 

REAL ZNEW 


NOW INTEGERS. 
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Si ee! Gn ep 1 a 


C 


C 


N IS THE MAIN LOOP GRIDPOINT COUNTER. 

INTEGER N 

COUNT AND COUNT1 ARE COUNTERS USED DURING AN AVERAGING PROCESS. 
INTEGER COUNT, COUNT1 

INDEX CORRELATES THE GRIDPOINT ALTITUDE WITH THE INDEX OF THE 
CORRECT DENSITY VALUE. 

INTEGER INDEX 

TICK COUNTS THE NUMBER OF OUTER LOOPS PERFORMED TO CONVERGENCE. 
INTEGER TICK 


FINALLY ARRAYS. 


DENSITY(S) IS THE DENSITY MEASURED AT 1000 FEET INTERVALS. 
REAL DENSTY(30) 

R(N) IS THE GRIDPOINT RADIAL POSITION IN FEET. 

REAL R(200) 

T(N) IS THE GRIDPOINT TENSION. THE ACTUAL GRIDPOINT LOCATION 
is AT N-1/2. 

REAL T(200) 

THETA(N) IS THE GRIDPOINT ANGULAR POSITION IN FEET. 

REAL THETA(200) 

Z(N) IS THE GRIDPOINT HEIGHT IN FEET. 

REAL Z(200) 


CEKKRARKKEKKKKKK KK KKK KKK KKK KK HHA KKK KK KK KKH KKH KKK HK KE KKK KEK KKK KEKE KEK KKK KK 


C 


SHAKE HANDS WITH THE OPERATOR. 
WRITE(6,*)! ' 

WRITE(6,*)! ' 

WRITE(6,*)! ' 

WRITE(6,*) ’ NAVAL POSTGRADUATE SCHOOL ' 
WRITE(6,*)' CLASS OF 1992 ' 
WRITE(6,*)’ ' 
WRITE(6,*)! \\ ' 
WRITE(6,*) ' Vac: > S ' 
WRITE(6,*)' * // ' 
WRITE(6,*) ' * ' 
WRITE(6,*) ' * ' 
WRITE(6,*) ' * ' 
WRITE(6,*)! * ' 
WRITE(6,*)! * ' 
WRITE(6,*)' * kK ' 
WRITE(6,*) ' é' 

WRITE(6,*)' TACL7A IS THE STEADY STEADY STATE SOLUTION TO THE’ 
WRITE(6,*)' ' TACAMO WIRE PROBLEM. ' 

WRITE(6,*)! ' 

PAUSE 


CRRA KAKAK KKK KKK KKK KK KKK KK KKK KH KK KK HK KK KK HHH KK KAK KKK AK KKK KK HK ke Kk Kk kek Kk Kk 


C 
C 


OPEN THE DATA FILES. 


OPEN (UNIT=11,FILE=’DATAO1.MAT’ ) 
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OPEN (UNIT=12, FILE='DATAO2 .MAT’ ) 
OPEN (UNIT=13, FILE='DATA03.MAT' ) 
OPEN (UNIT=14, FILE='DATA04.MAT’ ) 
OPEN (UNIT=15, FILE='DATA05.MAT’) 
OPEN (UNIT=16, FILE='DATA06.MAT'’ ) 
OPEN (UNIT=17, FILE='DATAO7.MAT’) 
OPEN (UNIT=18 , FILE='DATAO8 .MAT’ ) 
OPEN (UNIT=19, FILE='DATA09.MAT’ ) 
OPEN (UNIT=20, FILE='DATAO0O.MAT’) 
CEKK KKK KKK KK KKK KKK AKHKK KKK KKK KKK KKK HAHAHA HKKK HK HKKAK KKK HK KKK KKK HK KKK AK KKH 
FORMAT THE DATA FILES. 
DATAO].MAT IS UNFORMATTED. 


FORMAT (F12.6) 


FORMAT(F5.1) 
HAHAHAHAHAHA HA HA AHH HHA HRA HAHAH HHH HEH HEHEHE WDE EHH BEDE HEM A HHH HHH HHH HH 


INITIALIZE CONSTANTS. 


DROGUE CONSTANTS. 


63°) €) OO (NOOO 


AC=23.97/12.0 
6G=13,80/12.0 
CDD=0 .53 
CMACD=0 .0 
CLALDM=0. 61 
LEND=31.71/12.0 
SD=3.14159 
WD=81.95 


GENERAL CONSTANTS. 


‘o) 


G=32.174 

Pi=3. 1415926535879 
COUNT=1 

TICK=1 


WIRE CONSTANTS. 


OQ 


D=0.1582/12 
DELTAS=101.96 
MHU=0 .062107/G 


RR AND THTH ARE CONSTANTS USED TO START THE CALCULATIONS 
AT THE DROGUE. 


© 261.6) 


RR=1.0 
THTH=0 .001 


C READ THE DENSITY DATA FILE. 
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OPEN (31,FILE='’DENSTY.MAT’ ,STATUS='OLD’ , FORM=’FORMATTED’ , 


: ACCESS='SEQUENTIAL’ , RECL=5) 


READ (31,2,END=6)DENSTY 
GOTO 5 
CLOSE(31) 
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CALCULATE THE RADIUS OF THE TOWPLANE FLIGHT PATH AND THE THEDOT. 
THE PROGRAM WILL WORK AT ALL AIRSPEEDS AND ALTITUDES OF 
INTEREST AND FOR BANK ANGLES OF BETWEEN 5 AND 50 DEGREES. FOR 
BANK ANGLES ABOVE 45 DEGREES, AND AT TIMES AROUND THE LOCATION 
OF THE JUMP PHENOMENON, A SMALL ADJUSTMENT TO THE PROGRAM 

MAY BE REQUIRED. DOCUMENTATION WILL PRINT ON THE SCREEN IF 
THIS SITUATION OCCURS THAT WILL EXPLAIN HOW TO MAKE THE 
ADJUSTMENT. 

WRITE(6,*) INPUT AIRCRAFT KEAS, ALT IN FT, BANK ANGLE IN DEG’ 
READ(5,*) VEAS,ALTTP, PHI 

CONVERT KNOTS TO FEET PER SECOND. 

VEAS=VEAS*6076.1/3600.0 

CONVERT TO RADIANS. 

PHI=PHI*2 .0*P1I/360.0 

LOCAL GRIDPOINT DENSITY. A FILE OF MEASURED DATA MAY BE USED 
OR STANDARD ATMOSPHERE DATA MAY BE CALCULATED. COMMENT OUT THE 
METHOD NOT CHOSEN. 


INDEX=INT (ALTTP/1000.0)+1 
RHO=DENSTY ( INDEX) *0 .0023769/1013 .0 


RHO=0 .0023769%*(((518.69-0.0035662*ALTTP) /518 .69)** 
(Cie 0) * (1707 (-0.0035662%*53.3))+1.0)) 


RHO=0 .002378*(1-0.006875*Z(N) /1000.0)**4.256 
TRUE AIRSPEED FROM EQUIVALENT AIRSPEED. 
VTRUE=VEAS /SQRT (RHO /0 .0023769) 

CALCULATE TOWPLANE ORBIT RADIUS. 
RADTP=VTRUE** 2 /(G*SQRT(1/(COS (PHI) )**2-1)) 
ORBIT RATE CALCULATION. 


THEDOT=VTRUE/RADTP 
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THE INITIAL POSITION AT GRIDPOINT 1. 

THIS IS AN INITIAL GUESS AT THE BOTTOM POSITION. THIS VALUE IS 
ITERATED TO MATCH THE BOUNDARY CONDITION POSITION AT THE 
TOWPLANE. THESE VALUES MAY ALSO BE ADJUSTED TO FORCE THE 
MULTIPLE SOLUTIONS. TO FORCE MULTIPLE SOLUTIONS THE INITIAL 
POSITION MAY BE CHOSEN AT OR OUTSIDE THE TOWPLANE RADIUS OF 
ORBIT. 


THETA(1)=0.0 
IF (PHI .GT. 32.0*2.0*PI/360.0) THEN 
Z(1)=0.3*ALTTP 

Z(1)=0.4*ALTTP 
R(1)=0.2*RADTP 

R(1)=1.6*RADTP 

ELSE 

Z(1)=0,.6*ALTTP 

Z(1)=0.2*ALTTP 
R(1)=0.6*RADTP 

R(1)=1.2*RADTP 

ENDIF 
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CONTINUE 

AT THE DROGUE. 

GUESS AT THE POSITION USING THE FORCING FACTORS. 
R(2)=R(1)+RR 

THETA(2)=THETA(1)+THTH 
Z(2)=Z(1)+SQRT (DELTAS**2 -R(2)**2-R(1) *¥*2+ 
2*R(2)*R(1)*COS (THETA (2) - THETA(1) ) ) 


DENSITY AT GRIDPOINT 1. AGAIN, CHOOSE TABLE LOOKUP 
OR STANDARD ATMOSPHERE. 


INDEX=INT(Z(1)/1000.0)+1 
RHO=DENSTY ( INDEX) *0 .0023769/1013.0 


RHO=0 .0023769%*(((518.69-0.0035662*Z(1))/518.69) ** 
((-1.0)*(1.07(-0.0035662%5375))4.3)) 


RHO=0 .002378*(1-0.006875*Z(N)/1000.0)**4. 256 
LOCAL DYNAMIC PRESSURE. 


Q=0 . 5*RHO*(R(1)*THEDOT) **2 
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CALCULATE THE AOA AND SIDESLIP ANGLE OF THE DROGUE BY CALCU- 
LATING THE MOMENT EQUATION RESIDUES AND CHOOSING THE AOA AND 
SIDESLIP ANGLE AT WHICH THE RESIDUE COMES CLOSEST TO VANISHING. 
ONE DEGREE INCREMENTS ARE CONSIDERED SUFFICIENT. 


INITIALIZE THE RESIDUE HOLDERS SINCE THEY ARE REUSED AT EACH 
NEW BOUNDARY SHOT. 


STUFF2=1000 .0 
STUFF5=1000 .0 


LOOP THROUGH EACH ANGLE FROM O TO 90 DEGREES. 
NEGLECT ALL OTHER ANGLES THROUGH PHYSICAL REASONING. 


DO 50 N=1,90 
USE INDEX TO GET ANGLE IN RADIANS. 


ALFAD=N*2 .0*PI/360.0 
BETA=ALFAD 


CALCULATE THE LIFT CURVE SLOPE AT THIS ANGLE. 


CLALD=CLALDM 
NOW THE RESIDUES AT THIS ANGLE. 


FOR THE ANGLE OF ATTACK. 


FIRST CALCULATE THE PRODUCT OF CLALD AND ALFAD ACCOUNTING 
FOR THE EFFECTS OF STALL. 


STUFF 7=CLALD*ALFAD 

IF (ABS(ALFAD) .GT. 0.786)THEN 
STUFF7=0 .42* (ALFAD) / (ABS (ALFAD) ) 
ELSE 

ENDIF 

IF (ABS(ALFAD) .GT. 1.05) THEN 
STUFF7=0 . 2*(ALFAD) /(ABS(ALFAD) ) 
ELSE 

ENDIF 


STUFF 1=- STUFF7*Q*SD*AC*COS (ALFAD) - CDD*Q*SD*AC*SIN(ALFAD) + 


; WD*CG*COS (ALFAD ) +CMACD*Q*SD*LEND - CDD*SIN (BETA) *AC*Q*SD 


FOR THE SIDESLIP ANGLE. 


AGAIN CALCULATE THE PRODUCT OF CLALD AND BETA ACCOUNTING FOR 
THE EFFECTS OF STALL. 


STUFF /=CLALD*BETA 
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IF (ABS(BETA) .GT. 0./786)THEN 
STUFF/=0 .42* (BETA) / (ABS (BETA) ) 
ELSE 

ENDIF 

IF (ABS(BETA) .GT. 1.05) THEN 
STUFF/=0 .2* (BETA) / (ABS (BETA) ) 
BESt 

ENDIF 


STUFF4=WD*THEDOT**2*R (1) *CG*COS (BETA) /G- STUFF7*Q*SD*AC* 
COS (BETA)+CMACD*Q*SD 


SAVE RESIDUE AND ANGLE IF IT IS LESS THAN THE LAST. 
FOR THE ANGLE OF ATTACK. 


IF (ABS(STUFF1) .LT. ABS(STUFF2)) THEN 
STUFF2=STUFF1 

STUFF3=ALFAD 

ENDIF 


FOR THE SIDESLIP ANGLE. 


ELIMINATE THE POSSIBILITY OF MULTIPLE SOLUTIONS BEYOND 45 
DEGREES . 


IF (BETA .GT. PI/4.0) THEN 

GOTO 50 

ENDIF 

IF (ABS(STUFF4) .LT. ABS(STUFF5)) THEN 
STUFF5=STUFF4 

STUFF6=BETA 

ENDIF 

CONTINUE 


SELECT THE MINIMUM RESIDUE ANGLE OF ATTACK AS THE DESIRED VALUE. 


ALFAD=STUFF3 
BETA=STUFF6 


CALCULATE THE LIFT FORCE OF THE DROGUE DEFINED PERPENDICULAR TO 
THE RELATIVE FLOW AND CALCULATE THE LIFT AND SIDEFORCE OF THE 
DROGUE. 


CLALD=CLALDM* (PI -2.0*ALFAD) /PI 
LD=CLALD*ALFAD*Q*SD 
CLALD=CLALDM* ( PI-2.0*BETA) /PI 
ASFD=CLALD*BETA*Q*SD 


CALCULATE THE DRAG FORCE OF THE DROGUE DEFINED PARALLEL TO 
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THE RELATIVE FLOW. 


DD=CDD*Q*SD 
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NOW GET TENSION AT 1 AND POSITION AT 2 BY ASSUMING THE 
AERODYNMIC AND INERTIAL FORCES ON THE FIRST SEGMENT OF WIRE 
ARE SMALL COMPARED TO THE DROGUE FORCES. 


THE TENSION AT GRIDPOINT 1 MUST BALANCE THE VECTOR SUM OF THE 
FORCES UPON THE DROGUE IN ORDER TO SATISFY STATIC EQUILIBRIUM. 


STUFF1=WD*THEDOT**2*R (1) /G 
T(1)=SQRT( (STUFF1 - ASFD) **2+DD**2+(LD-WD)**2) 


ASSUME THAT THE FORCES ON THE FIRST SEGMENT OF THE WIRE ARE 
SMALL COMPARED TO THE FORCES ON THE DROGUE. 


T(2)=T(1) 


THE RESULTANT FORCE MUST ALSO BE IN THE RECIPROCAL DIRECTION 
TO MAINTAIN STATIC EQUILIBRIUM. 


R(2)=R(1) - (STUFF1-ASFD)*DELTAS/T(1) 
THETA (2 )=THETA(1)+DD*DELTAS /(T(1)*((R(1)4+R(2))/2)) 


USE THE TWO POINT SEGMENT LENGTH CONSTRAINT TO GET THE Z AT 
GRIDPOINT 2. 


Z(2)=Z(1)+SQRT (DELTAS **2-R(2) **2-R( 1) **2+ 
2*R(2)*R(1) *COS (THETA (2 ) - THETA(1) ) ) 
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ITERATE TO THE TOP. 

MAKE INITIAL GUESSES AT THE ITERATION VALUES FOR GRIDPOINT TWO. 
Al=T(1)*(R(2)-R(1))/DELTAS 
Bl=T(1)*(R(2)+R(1) )*(THETA(2) -THETA(1) ) /(2*DELTAS ) 
Cl=T(1)*(Z(2)-Z(1))/DELTAS 

DO 3500 N=3,200 


GUESS AT NEXT POINT. USE THE THREE POINT UNIT TANGENT VECTOR 
CONSTRAINT TO GET THE Z AT GRIDPOINT N. 


R(N)=R(N-1)+R(N-1) -RCN- 2) 

THETA (N)=THETA(N- 1)+THETA(N-1) -THETA(N-2) 
Z(N) =4*DELTAS**2- (R(N) -R(N- 2) )**2-R(N-1)**2* 
(THETA (N) - THETA(N-2) )**2 
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PROTECT AGAINST A NEGATIVE RADICAND DURING THE ITERATIONS. 


IF (Z(N) .LT. 0.0) THEN 
Z(N)=0.000001 

ELSE 

ENDIF 
Z(N)=Z(N-2)+SQRT(Z(N) ) 


T(N)=T(N-1)+MHU*DELTAS*G/(R(N-1)*(THETA(N-1)-THETA(N-2)) /DELTAS ) 
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THE INNER LOOP FOR MOVING FROM ONE GRID POINT TO THE NEXT. 
INITIALIZE COUNTERS FOR LATER USE. 


COUNT=1 
COUNT1=1 


INITIALIZE SUMMERS FOR USE IN AVERAGING LATER. 
A2A2=0.0 

B2B2=0.0 

C2C2=0 .0 

CONTINUE 

CALCULATE DENSITY. 


INDEX=INT(Z(N) /1000.0)+41 
RHO=DENSTY ( INDEX) *0 . 0023769/1013.0 


RHO=0 .0023769*(((518.69-0.0035662*Z(N) )/518.69)** 
((-1.0)*(1.0/(-0.0035662*53.3))+1.0)) 


RHO=0 .002378%*(1-0.0068/5*Z(N) /1000.0)**4 .256 

INCREMENT COUNTER. 

COUNT=COUNT+1 

NOW THE ITERATIONS. 

CALCULATE SOME TRANSCENDENTALS USED IN LATER CALCULATIONS. 
STUFF1=R(N-1)*(THETA(N) -THETA(N- 2) )/(2*DELTAS ) 
STUFF2=(1-STUFF1**2) 


STUFF3=3*SQRT (STUFF2 ) -4* (SQRT (STUFF2 ) ) **3 


CALCULATE CN, CA AND CY. 


108 


C3°C? ©) ©) 


CN=1 .09*STUFF2-0.08*4 .O*STUFF2*STUFF1**2 
TPS CotUrrieer7Lr, OF7e7)” THEN 
CA=0 .0 

ELSE 

CA=0.017 

ENDIF 

PreecstUPRI bb. 075) THEN 
CY=0.0 

ELSE 

CY=0.0148*STUFF3*3 

ENDIF 

CY=0 .0 


THE ITERATIVE EQUATIONS. 


STUFF1=(T(N)+T(N-1))*R(N-1)*(THETA(N) - THETA(N-2) )**2 


: /(8*DELTAS**?2) 


STUFF2=CN*R(N-1)*(THETA(N) - THETA(N-2))*(R(N) -R(N-2))/ 
(4*DELTAS**2) 
STUFF3=CA*R(N-1)*(THETA(N) -THETA(N-2) )*(R(N)-R(N-2))/ 
(4*DELTAS**2) 

STUFF5=(R(N)-R(N-2) )**2+R(N-1)**2* (THETA (N) - THETA(N-2) ) **2 
STUFF5=SQRT(STUFF5) 

STUFF4=R (N-1)*(THETA(N) - THETA(N-2) )*0. 5*RHO*(R(N-1) *THEDOT) **2 


: *D*CY/STUFF5 


A2=A1+(STUFF1-0.5*RHO*(R(N-1)*THEDOT) **2*D* (STUFF2 - STUFF3 ) 
- MHU* THEDOT**2*R (N- 1)+STUFF4 ) *DELTAS 


STUFF1=(T(N)+T(N-1))*(R(N) -R(N-2) )*(THETA(N) - THETA(N- 2) ) 


: /(8*DELTAS**2) 


STUFF2=CN* ((R(N-1)**2%* (THETA (N) -THETA(N-2) ) **2/ 
(4*DELTAS**2) )-1) 
STUFF3=CA*R(N-1)*(THETA(N) - THETA(N-2) ) /(2*DELTAS) 
STUFF4=(R(N) -R(N-1))*0.5*RHO*(R(N-1)*THEDOT )**2*D*CY /STUFF5 
B2=Bl - (STUFF1+0 . 5*RHO*(R(N- 1) *THEDOT) **2*D* (STUFF2 - STUFF3) 


: +STUFF4)*DELTAS 


STUFF1=CN*R(N-1)*(THETA(N) - THETA(N-2))*(Z(N) -Z(N-2)) 


: /(4*DELTAS**2 ) 


STUFF2=CA* (Z(N) -Z(N-2) ) /(2*DELTAS) 
C2=Cl+(-0.5*RHO*(R(N-1)*THEDOT) **2*D* (STUFF1 - STUFF2 ) 


>: +MHU*G)*DELTAS 


QO 


T(N)=SQRT (A2**24+B2**2+C2**2 ) 
R(N)=R(N-1)+A2*DELTAS/T(N) 
THETA(N)=THETA(N-1)+B2*2*DELTAS /(T(N)*(R(N)+R(N-1))) 
Z(N)=Z(N-1)+C2*DELTAS /T(N) 


HERE, AGAIN APPLY A LITTLE TRICK. FOR 8088 SYSTEMS AND 
ON VERY RARE OCASSIONS FOR 16 BIT SYSTEMS, THE PROGRAM TENDS TO 
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ENTER A LIMIT CYCLE LIKE BEHAVIOR IN TENSION AFTER 100-200 
ITERATIONS WHICH CYCLES EVERY 2-3 STEPS. 

IT IS SOLVED BY SIMPLY AVERAGING THE VALUES OVER 20 

STEPS ANY TIME 400 ITERATIONS ARE EXCEEDED. THIS SIMPLE 
SOLUTION HAS BEEN FOUND TO WORK FOR ALL CASES CHECKED. 
AGAIN, USING THE 16 BIT 486 TYPE MACHINE TO RUN THE SOFTWARE 
PRETTY MUCH PRECLUDES THE PROBLEM EXCEPT IN THE RAREST OF 
CIRCUMSTANCES. 
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IF (COUNT .LT. 400) THEN 
GOTO 1200 
ELSEIF (COUNT1 .LT. 21) THEN 
COUNT1=COUNT1+1 
A2A2=A2A2+A2 
B2B2=B2B2+B2 
C202=62624+@2 
GO TO 1200 
ELSE 
A2=A2A2/20.0 
B2=B2B2/20.0 
C2=C2G2 2000 
GOTO 3000 
ENDIF 
1200 CONTINUE 


C COMPARE TO THE CONVERGENCE CRITERIA. 


IF (ABS((A2-A22)/A2) .GT. 0.0001) THEN 
GOTO 2000 

ELSEIF (ABS((B2-B22)/B2) .GT. 0.0001) THEN 
GOTO 2000 

ELSEIF (ABS((C2-C22)/C2) .GT. 0.0001)THEN 
GOTO 2000 

ELSE 

GOTO 3000 

ENDIF 


2000 CONTINUE 
A22=A2 
B22=B2 
C22=C2 
GOTO 1000 


3000 CONTINUE 
C MOVE THE ITERATION VALUES FORWARD ONE TIME STEP. 
C 

Al=A2 


B1l=B2 
C1l=C2 
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UPDATE THE POSITION VALUES AND THE TENSION USING ALL THE 
LATEST INFORMATION. 


T(N)=SQRT (A2**24+B2**24C2¥**2 ) 
R(N)=R(N-1)+A2*DELTAS /T(N) 

THETA (N)=THETA (N-1)+B2*2*DELTAS /(T(N)*(R(N)+R(N-1))) 
Z(N)=Z(N-1)+C2*DELTAS /T(N) 


RESET THE SUMMERS FOR THE NEXT LOOP. 


A2A2=0.0 
B2B2=0.0 
C2C2=0.0 
CONTINUE 
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MATCH THE BOUNDARY CONDITION AT THE TOP. THAT IS, MATCH THE 
TOP POSITION OF THE WIRE TO THAT OF THE TAIL OF THE TOWPLANE. 
THE ERROR AT THE TOP IS INTERPOLATED AND USED AS A FACTOR TO 
ADJUST THE BOTTOM POINT. ON OCCASSION, THE MULTIPLICATIVE 
FACTORS BELOW MUST BE ADJUSTED TO ENSURE CONVERGENCE OF THE 
BOUNDARY CONDITION SHOOTING ROUTINE. WHEN THIS IS REQUIRED, 

THE PROGRAM PROVIDES GUIDANCE ON HOW TO PERFORM THE ADJUSTMENTS. 


MATCH THE Z AT EACH CHOSEN RADIUS. 
IF (ABS(Z(200)-ALTTP) .GT. 50.0) THEN 


IF (R(1) .LT. 700.0) THEN 
ZNEW=Z(1)+0.1* (ALTTP-Z(200) ) 
GOTO 3550 

ELSEIF (R(1) .LT. 1000.0) THEN 
ZNEW=Z(1)+0.2*(ALTTP-Z(200) ) 
GOTO 3550 

ELSEIF (R(1) .LT. 1700.0) THEN 
ZNEW=Z(1)+0.25*(ALTTP-Z(200) ) 
GOTO 3550 

ELSE 
ZNEW=Z(1)+0.3*(ALTTP-Z(200) ) 
ENDIF 


CONTINUE 
NOW ADJUST TO A NEW RADIUS. 
Z(1)=ZNEW 


RNEW=R (1) 
GOTO 3600 
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3650 


ELSEIF (ABS(R(200)-RADTP) .GT. 10.0) THEN 


IF (R(1) .LT. 800.0) THEN 
RNEW=R(1)+0.07* (RADTP-R(200) ) 
GOTO 3560 

ELSEIF (R(1) .LT. 1000.0) THEN 
RNEW=R(1)+0.15*(RADTP-R(200) ) 
GOTO 3560 

ELSEIF (R(1) .LT. 2000.0) THEN 
RNEW=R(1)+0.15*(RADTP-R(200) ) 
GOTO 3560 

ELSE 

RNEW=R(1)+0.3* (RADTP-R(200) ) 
ENDIF 


CONTINUE 


R(1)=RNEW 
ZNEW=Z (1) 
GOTO 3600 
ELSE 

GOTO 3700 
ENDIF 


CONTINUE 


R(1)=RNEW 
Z(1)=ZNEW 


CONTINUE 


PRINT SOME OUTPUTS USEFUL IN KEEPING UP WITH HOW THE BOUNDARY 
CONDITION SHOOTING ITERATIONS ARE PROGRESSING. 


WRITE(6,*) °RCL),ZC1)=" Ri) Zo 

WRITE(6,*) ‘R(200) ,Z(200) ,T(200)=' ,R(200) ,Z(200) , T( 200) 
WRITE(6,*) ‘RADTP,ALTTP=’ ,RADTP,ALTTP 

WRITE (6, ¥) 1 8 KKK HAA HHH HAA HAHAHA HE HAHA HHA HEE HA EER EEE EEE TEE IEE TEE IEE TE! 


JUST BOOKKEEPING HERE FOR THE NEXT LOOP. 


DO 3650 N=2,200 
R(N)=0.0 
THETA(N)=0.0 
Z(N)=0.0 
T(N)=0.0 


CONTINUE 
TICK=TICK+1 
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CHECK TO SEE IF THE UPDATE MULTIPLIERS ABOVE HAVE PLACED THE 
BOUNDARY SHOOTING PROCESS INTO AN "INFINITE LOOP” AND IF SO 
STOP THE PROCESS AND TELL THE OPERATOR HOW TO FIX THE PROBLEM. 


IF (TICK .GT. 400) THEN 
WRITE(6,*)'YOU HAVE EXCEEDED 400 ATTEMPTS TO SHOOT THE BOUNDARY’ 
WRITE(6,*)'CONDITION AT THE TOP. TO MAKE THE PROGRAM RUN ' 
WRITE(6,*)'CORRECTLY YOU MUST ADJUST THE MULTIPLIER’ 
WRITE(6,*)'COEFFICIENTS IN THE PROGRAM BETWEEN LINES 3500 AND’ 
WRITE(6,*)'3560. NOTE THE PRINTOUT OF R(1) AND Z(1) ABOVE AS’ 
WRITE(6,*)'THEY FLOP BACK AND FORTH ON EITHER SIDE OF THE ' 
WRITE(6,*) ‘CORRECT VALUE. FIRST TRY TO HALVE THE COEFFICIENT’ 
WRITE(6,*) 'ASSOCIATED WITH THE R(1) VALUES YOU SEE (NOTE THE’ 
WRITE(6,*)'R(1) .LT. STATEMENT). IF YOU ARE STILL NOT ' 
WRITE(6,*)'CONVERGING, TRY THE Z(1) COEFFICIENT. YOU WILL HAVE’ 
WRITE(6,*)'TO RE-COMPILE AFTER EACH FIX. I RECOMMEND THAT’ 
WRITE(6,*)’YOU WRITE DOWN THE ORIGINAL VALUES AND RETURN THEM’ 
WRITE(6,*)’WHEN YOU ARE DONE. THESE MULTIPLIERS GIVE YOU THE’ 
WRITE(6,*) ‘QUICKEST CONVERGENCE FOR THE MAJORITY OF CASES OF' 
WRITE(6,*)’ INTEREST. ’ 


GOTO 6000 
ELSE 
ENDIF 


GOTO 20 
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3700 


CONTINUE 
WRITE(6,*) ‘NUMBER OF OUTER LOOPS=' , TICK 
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WRITE OUTPUT TO THE THREE DATA FILES. 
FIRST, CONVERT SOME BACK TO MORE RECOGNIZABLE FORMS. 


VEAS=VEAS*3600.0/6076.0 
PHI=PHI*360.0/(2.0*PI) 


NOW, WRITE SOME OF THE CONSTANTS TO FILES. 


WRITE(11,*) ‘TOWPLANE KEAS=' ,VEAS, ‘KTS’ 

WRITE(11,*) ‘TOWPLANE ALTITUDE=' ,ALTTP, ' FEET’ 
WRITE(11,*) ‘TOWPLANE BANK ANGLE=' , PHI, ‘DEGREES’ 
WRITE(11,*) ‘TOWPLANE ORBIT RADIUS=' ,RADTP, ‘FEET’ 
WRITE(11,*) ‘TOWPLANE ORBIT RATE=' , THEDOT, 'RAD/SEC’ 
WRITE(11,*) ‘DROGUE ANGLE OF ATTACK=’ ,ALFAD, ‘RAD’ 
WRITE(11,*) ‘DROGUE SIDESLIP ANGLE=' ,BETA, 'RAD’ 
WRITE(11,*) ‘DROGUE ORBIT RADIUS=' ,R(1), ‘FEET’ 
WRITE(11,*) ‘DROGUE AIRSPEED=' , THEDOT*R(1) , ‘FT/SEC’ 
WRITE(11,*) '‘VERTICALITY=' , (ALTTP-Z(1))/(200.0*DELTAS) 
WRITE(11,*) ‘DROGUE/TOWPLANE SEPARATION=' ,ALTTP-Z(1), ‘FEET’ 
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4000 


STUFF1=ASIN((R( 199) -R( 200) ) /DELTAS ) 
STUFF1=STUFF1*57. 3 

WRITE(11,*) ‘WIRE TRAIL ANGLE=', STUFF1 
WRITE(16,1) THEDOT 

WRITE(17,1) CLALD*ALFAD 

WRITE(19,1) VTRUE 

WRITE(20,1) PHI*2.0*PI/360.0 


WRITE THE POSITION AND TENSION TO DATA FILES. 


DO 4000 I=1, 200 
WRETE( 1271) 3G 
WRITE(13,1) THEDACI) 
WRITE(14,1) Z(I) 
WRETEC ISS Deo 
CONTINUE 
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C 


C 
C 
C 


5000 


BUILD A FILE THAT CONTAINS THE ANGLE OF ATTACK OF EACH WIRE 
GRIDPOINT. NOTE THAT THIS DOES NOT INCUDE THE FIRST OR LAST 
POINTS. 


DO 5000 I=2,199 

STUFF1=ACOS (R(1I)*(THETA(I+1) - THETA(I-1)) /(2*DELTAS) ) 
STUFF1=STUFF1*360/(2*PI) 

WRITE(18 ,1)STUFF1 

CONTINUE 
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C 
C 
C 
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5100 


THESE FILES CONTAIN THE REYNOLDS NUMBER AND TRUE AIRSPEED AT 
EACH GRIDPOINT. 


OPEN (UNIT=70, FILE='VTRUE.MAT’ ) 

OPEN (UNIT=71,FILE='RE.MAT’ ) 

DO 5100 I=1,200 

INDEX=INT(Z(I) /1000.0)+1 

RHO=DENSTY ( INDEX) *0 .0023769/1013.0 

RHO=0 .002378*(1-0.006875*Z(N) /1000.0)**4. 256 
RHO=0 .0023769*(((518.69-0.0035662*Z(N) ) /518. 69) ** 
((-1.0)*(1.0/(-0.0035662*53.3))+1.0)) 

VTRUE=R(1)*THEDOT 

STUFF1=RHO*VTRUE*D/ .00000038 

WRITE(70,1)VTRUE 

WRITE(71,1)STUFF1 

CONTINUE 

CLOSE(70) 

CLOSE(71) 
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WRITE(6,*)! ' 

WRITE(6,*) ‘RUN COMPLETE!!!’ 

WRITE(6,*)’ ' 

WRITE(6,*) ‘THE OUTPUT IS INCLUDED IN THE FOLLOWING: ' 
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WRITE(6,*)’DATAO1. 
WRITE(6,*) ’DATAO2. 
.MAT CONTAINS EACH GRIDPOINT THETA COORD. ' 
MAT CONTAINS EACH GRIDPOINT Z COORD. ’ 


WRITE(6,*) ‘DATAO3 
WRITE(6,*) 'DATAO4 


WRITE(6,*) 'DATAOS. 
.MAT, DATAO7.MAT AND DATAOO.MAT CONTAIN VALUES’ 


WRITE(6,*) ’DATAO6 


MAT CONTAINS MISCELLANOUS VALUES OF INTEREST. ’ 
MAT CONTAINS EACH GRIDPOINT RADIAL COORD. ’ 


MAT CONTAINS EACH GRIDPOINT TENSION VALUE. ’ 


WRITE(6,*) ’REQUIRED BY THE DYNAMIC SOLUTION PROGRAM. ’ 


WRITE(6,*) 'DATAO8. 


MAT CONTAINS THE TRUE ANGLE OF ATTACK OF’ 


WRITE(6,*) ’GRIDPOINTS 2 THROUGH 199. VTRUE .MAT CONTAINS THE’ 
WRITE(6,*)’TRUE AIRSPEED AT EACH GRIDPOINT AND RE.MAT CONTAINS’ 


WRITE(6,*) ’THE REYNOLDS NUMBER AT EACH GRIDPOINT. ’ 
CREE EE HEHEHE HEHE DE DE HA HAR WE HEA IEEE WHEE EEE BER HE EEE BEE EEE EE EE ETE ETE DE DEE TEBE IE 


6000 CONTINUE 
END 
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